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Preface 


The  purpose  of  this  st«Hy  was  to  investigate  the  effect 
of  using  a  Poisson  approxinatinn  in  the  WARS  ASSESS  subrou¬ 
tine.  This  approximation  is  a  snail  part  of  the  ASSESS  sub¬ 
routine  which  in  turn  is  a  small  part  of  the  WARS  model. 
Our  altn  was  not  to  evaluate  the  preformance  of  WARS  or 
ASSESS  but  only  to  study  the  effect  of  this  particular 
approximation.  We  hope  that  in  some  small  way  we  have  added 
to  the  understanding  of  this  important  model. 


V’e  would  like  to  thank  our  advisors,  Mai  Joseph  W. 
Coleman  and  Lt  Col  James  N.  Bexfield  of  the  Air  Force  Insti¬ 
tute  of  Technology  at  Wr igh t-Pa t t e r son  AFB,  Ohio,  who  have 
guided  this  study  to  completion.  Sincere  gratitude  is  also 
expressed  to  Mai  Pouglas  Rippy  and  Capt  Steven  Schroeder  for 
their  expert  council.  Finally  we  would  like  to  acknowledge 
the  patience  of  our  wives  and  families  in  this  demanding 
undertak ing . 

Steven  W.  Weiss 
Neil  A.  Youngman 
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Abstract 

The  purpose  of  this  study  was  to  determine  the  relative 
effect  of  using  the  Poisson  approximation  in  the  Wartime 
Assessment  and  Requirements  System  (WARS)  assessment  subrou¬ 
tine  and  to  gain  a  better  understanding  of  the  exact  error- 
causing  mechanisms  involved. 

To  accomplish  this  two  Fortran  computer  programs  were 
developed,  one  to  compute  the  expected  number  of  not  mission 
capable  aircraft  using  the  accurate  mathematical  calcula¬ 
tions,  the  other  to  compute  the  same  figure  using  the  Pois¬ 
son  approximation.  These  programs  were  used  to  evaluate  the 
approximation  caused  error  for  different  parts  hierarchies 
and  data  sets  . 

The  analysis  identified  two  distinct  causes  for  the 
error  induced  by  the  approximation.  These  sources  of  error 
were  confirmed  by  running  test  cases  specifically  tailored 
to  eliminate  the  error-causing  characteristics  and  noting 
that  no  approximation  error  resulted. 

The  approximation  error  was  found  to  fluctuate  in  sign 
and  magnitude  for  different  cases.  Sensitivity  analysis  was 
performed  to  identify  the  sensitive  parameters. 


THE  EFFECT  OF  THE 


POISSON  APPROXIMATION  ON  THE 
WARTIME  ASSESSMENT  AND  REQU  IP  EMF,  NTS  SYSTEM 
ASSESSMENT  RESULTS 

I  Introduction 


Background 

The  War  Reserve  Materiel  Concept .  In  peacetine  an 
operational  unit  receives  its  required  spare  parts  through 
normal  supply  channels;  however,  during  contingencies  the 
unit  may  be  deployed  to  another  location  and  be  unable  to 
maintain  its  established  supply  lines*  If  the  nohilization 
is  in  response  to  a  contingency  plan,  logistical  support  is 
provided  through  the  War  Reserve  Materiel  (WRM)  program, 
which  is  designed  to  support  the  wartime  mission  until  pro¬ 
duction  and  airlift  can  ensure  adequate  resupply.  This  ini¬ 
tial  resupply  is  accomplished  by  stockpiling  war  materiel 
and  prepositioning  a  portion  of  it  in  the  theatre  of 
intended  use  prior  to  the  beeinning  of  hostilities  (Ref  1: 
Ch  1,  1). 

WRM  is  divided  into  three  separate  areas.  War  Readiness 
Spares  Kits  (WRSK),  Base  Level  Self-sufficiency  Spares 
(BLSS),  and  Other  War  Reserve  Materiel  (OWRM).  WRSK  is  a 


special  category  of  WRM  which  Is  defined  as: 


An  air  transportaM  r  r.irlMpc  of  V  ;•«  ”  spares,  repair 
parts  and  related  maintenance  supplies  required  to  sup¬ 
port  planned  wartime  or  contingency  operations  of  a 
weapon  or  support  system  for  a  specifier?  period  of  time 
pending  resupply  (Ref  2:  Ch  14,  1CA). 

Current  Air  Force  directives  normally  specify  a  30-dav  time 
period  without  resupply.  WRSKs  are  authorized  for  specific 
units  which  have  weapon,  support,  or  mobile  command  and  con¬ 
trol  systems,  and  are  normally  prepositioned  with  the  using 
unit  (Ref  3:  Ch  2,  2).  Although  WRM  consists  of  various 
types  of  materiel  such  as  petroleum,  munitions,  rations, 
spare  parts,  and  communication  eouipment,  WRSKs  contain  only 
mission  essential  replacement  parts  (Ref  2:  Ch  14,  10A). 
BLSS  is  similar  to  WRSK  in  that  its  purpose  is  also  to 
insure  an  adequate  supply  of  spare  parts  for  units  during 
wartime  operations,  however,  unlike  WRSK,  "BLSS  is  a  non- 
mobile  spares  package  that  augments  inplace  operating 
spares"  (Ref  3s  Ch  2,  1)  at  selected  forward  operating 
bases.  OWRM  is  WRM  that  is  not  included  in  WRSK  or  BLSS. 

Current  Requirements  Methodology.  The  basic  WRM  con¬ 
cept  is  easily  understood;  however,  the  actual  determination 
of  "which  parts"  and  "how  many”  to  include  in  a  WRM  package 
is  a  difficult  task.  The  large  number  of  variables  and  the 
uncertainties  associated  with  them  add  to  the  complexity  of 
the  problem.  Determining  optimal  quantities  is  extremely 
important  because  in  time  of  war,  shortages  of  essential 


spare  parts  can  severely  degrade  mission  accomplishment. 
Conversely,  stockpiling  excess  quantities  or  unnecessary 


parts  is  very  expensive,  especially  considering  the  fact 
that  the  current  WRSK  inventory  alone  is  valued  at  near  S3. 6 
billion  (Ref  4).  One  can  easily  see  that  ninor  inprovements 
in  methodology  can  result  in  na.jor  improvements  in  nission 
capability  and  substantial  dollar  savings. 

WRSK,  BLSS,  and  OWRM  quantities  are  computed  by  similar 
methods,  so  for  the  sake  of  simplicity  the  subsequent  dis¬ 
cussion  will  be  limited  to  WRSK. 

Currently  the  kits  are  built  using  a  "WRSK  candidate 
list.”  This  is  a  unique  listing  for  each  aircraft  type  or 
other  system,  which  has  been  developed  over  the  years  by 
users  and  system  managers  to  specify  which  spare  parts  are 
mission  essential  and  are  therefore  authorized  for  WRSK. 
This  compilation  is  then  used  as  a  "shopping  list"  of  possi¬ 
ble  items  for  inclusion  in  the  kit. 

To  answer  the  "how  many  of  each  item"  question,  the  Air 
Force  Logistics  Command  has  developed  several  different 
analytical  techniques.  Two  techniques  currently  used  are 
referred  to  as  the  conventional  method  and  the  marginal 
analysis  method  and  are  described  below. 

A  WRSK  built  by  the  conventional  method  (called  a  con¬ 
ventional  kit)  basically  includes  quantities  equal  to  the 
maximum  expected  shortage  for  each  spare  part  for  the  entire 
WRSK  period.  These  Quantities  are  determined  in  the  follow- 


1  np  manner.  I’sine  historical  data,  a  mean  failure  rate  per 
flying  hour  is  calculated  for  each  item  by  dividing  the 
total  number  of  failures  in  a  certain  time  period  by  the 
number  of  flyinp  hours  in  the  period.  This  mean  failure 
rate  is  multiplied  by  the  number  of  items  per  aircraft  and 
the  expected  wartime  Daily  Flyinp  Hour  Propram  (DFHP)  to 
arrive  at  an  expected  number  of  failures  for  each  day.  If 
no  ma intenance  capability  exists  or  if  the  itc  is  a  non¬ 
reparable  type,  the  maximum  expected  shortage  for  he  period 
is  simply  the  sum  of  the  daily  expected  failures  owever, 
if  maintenance  capability  does  exist,  the  number  of  repar¬ 
able  items  repaired  on  each  day  is  computed  using  an  average 
repair  rate  (called  Base  Repair  Cycle  or  BRC)  and  an  initial 
setup  time  for  a  maintenance  shop.  Knowing  the  expected 
number  of  failures  for  each  day  and  the  expected  number  of 
repaired  items  each  day,  the  maximum  expected  shortage  can 
he  found.  This  maximum  expected  shortage  is  the  quantity  of 
that  particular  Item  included  in  the  conventional  kit.  The 
process  is  then  repeated  for  each  item  on  the  WRSK  candidate 
list  for  the  weapon  system  under  consideration  (Ref  5:  Ch  2, 
1-2;  6:  1-8). 

The  Marginal  Analysis  (MA)  method  is  designed  to  create 
a  new  WRSK,  with  the  "same  performance"  as  the  conventional 
WRSK,  at  a  cheaper  price.  Performance  is  defined  in  terms 
of  two  parameters,  maximum  expected  WRSK  shortage  and  max¬ 
imum  expected  number  of  Not  Mission  Capable  (NMC)  aircraft, 


both  taken  over  the  desired  number  of  davs.  for  example,  If 
the  URSK  neriod  to  be  supported  Is  15  days  and  a  tertian  kit 
is  evaluated  as  Indicated  In  the  follow  in  p  table,  then  the 
maximum  expected  number  of  MMO  aircraft  is  3.1  and  the  max¬ 
imum  expected  part  shortage  is  4.2.  These  two  figures 
represent  the  performance  of  this  particular  kit. 


EVALUATION  OF  WRSK  X 

DAY  EXPECTED  MO.  OF  PARTS  SHORT  EXPECTED  MO.  OF  A/C  NMC 
1  0  0 

2  0  0 

3  0  .7 

4  1.1  1 

5  2  1.3 

6  2.3  1.5 

7  3  2.1 

8  3.2  2.4 

9  4  2.9 

10  4.2  3.1 

11  3.8  2.3 

12  3.3  2.1 

13  3  1.9 

14  2.4  1.5 

15  2  1.4 

Table  1-1 

Note:  The  decrease  after  day  10  could  he  caused  bv  repair 

capability  becoming  operational  or  a  decrease  in  the  DFHP. 


We  will  use  aircraft  for  illustrative  purposes  although 
WRSK  is  not  specifically  limited  to  aircraft  and  flying 
units  . 

An  MA  WRSK  is  determined  by  first  building  a  conven- 


6 


tional  kit  and  evaluating  it  with  respect  to  the  shortage 


and  KMC  parameters,  that  is,  determining  the  maximum 
expected  number  of  broken  aircraft  and  the  maximum  expected 
number  of  parts  short,  over  the  period  for  this  kit.  These 
parameters  are  then  used  as  lower  bounds  for  the  MA  kit, 
which  is  huilt  from  a  small  initial  base,  normally  adding 
items  one  at  a  time.  The  next  item  to  be  added  is  the  one 
that  provides  the  greatest  reduction  in  the  parameters  per 
dollar  spent.  The  reduction  in  expected  shortage  due  to  the 
addition  of  a  given  part  to  a  WRSK  is  found  by  assuming  that 
failures  occur  according  to  a  Poisson  probability  distribu¬ 
tion  and  using  expected  value  calculations.  The  reduction 
in  expected  KMC  (EfNMC))  aircraft  is  found  in  the  same 
manner  after  making  the  assumption  of  full  cannibalization 
of  parts.  This  assumption  concentrates  all  the  shortages  or 
failures  into  the  minimum  number  of  aircraft  and  hence  sim¬ 
plifies  the  problem.  Parts  which  cannot  be  cannibalized  are 
not  considered  in  MA,  their  quantities  determined  solely  by 
the  basic  conventional  method.  This  process  of  adding  the 
most  cost  efficient  item  to  the  kit  is  continued  until  the 
kit  has  the  same  performance  as  the  conventional  kit  but  at 
less  cost  (Ref  5:  Ch  3,  1;  6). 

The  current  methodology  for  determining  WRSK  require¬ 
ments  has  three  inherent  problem  areas.  The  first  shortcom¬ 
ing  is  the  treatment  of  cannibalization.  In  the  present 
methodology  full  cannibalization  is  assumed  and  anv  part 
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which  is  cons  id  p  red  "non-cannahlp  is  not  i  nc  1  ml  p<i  in  the 
optimization  process.  The  WRSK  quantities  for  these  non- 
cannable  parts  are  determined  by  the  conventional  method. 
This  procedure  results  in  the  non-cannable  portion  of  a  WRSK 
being  sub-optimal  and  therefore  more  costly  than  needed  (Ref 

7). 

A  second  problem  is  the  relationship  between  Line 
Replaceable  Units  (LRUs)  and  Shop  Replaceable  Units  (SRUs). 
An  SRU  is  a  component  part  of  an  LRU  but  because  of  its 
inaccessibility  it  can  be  repaired  only  in  the  maintenance 
shop.  For  example,  the  LRU  is  the  main  assembly  or  black 
box  which  can  be  replaced  on  the  flight  line  while  the  SRU 
is  the  subassembly  or  component  in  the  black  box  which  can 
be  replaced  only  in  the  shop.  Currently  SRUs  and  LRUs  are 
considered  as  unrelated  parts  whereas  they  should  have  some 
sort  of  indenture.  relationship.  An  SRU  should  not  fail 
without  its  corresponding  LRU  also  failing.  For  example,  if 
an  amplifier  in  a  radio  fails,  the  entire  radio  should  also 
fail.  The  present  methodology  does  not  take  this  into 
account.  The  fact  that  LRUs  are  themselves  a  source  of 
parts  for  SRUs  is  also  not  addressed  (Ref  7). 

The  final  problem  area  is  inclusion  of  a  goal  which 
limits  the  number  of  backordered  parts  (called  stock  due  out 
or  SDO).  The  MA  method  builds  a  new  WRSK  with  the  same 
level  of  performance  as  the  conventional  kit  but  at  a 


cheaper  price  (Ref  5).  The  rationale  behind  this  method  is 
that  the  same  number  of  aircraft  will  he  N’MC  because  of  part 
shortages,  but  the  parts  short  will  be  expensive  parts 
rather  than  cheap  parts.  Performance  is  defined  in  terms  of 
two  parameters,  maximum  expected  SDO  and  maximum  expected 
number  of  NMC  aircraft,  both  taken  over  the  VRSK  period 
under  consideration.  The  requirement  that  the  M A  kit  not 
exceed  the  conventional  kit  SDO  goal  is  inappropriate 
because  an  SDO  goal  is  not  an  accurate  measure  of  a  VJHSK's 
ability  to  support  the  mission.  Imposing  a  maximum  SDO 
level  unduly  restricts  the  VRSK  and  costs  extra  dollars  (Ref 

7). 

The  Air  Force  realized  the  shortcomings  in  the  present 
methodology  and  began  development  of  the  Wartime  Assessment 
and  Requirements  System  (WARS)  to  improve  URSK  and  overall 
WRM  performance. 

Wartime  Assessment  and  Requirements  System  (Ref  fl )  . 
WARS  is  a  mathematical  optimization  model  used  to  determine 
the  quantities  of  WRM  items  in  a  kit  which  will  yield  the 
best  performance  per  dollar  spent.  The  Air  Force  Logistics 
Command  is  currently  developing  this  model  which  will  use  a 
MA  method  of  optimization  similar  to  the  method  previously 
described.  The  main  differences  are  that  WARS  uses  an  NMC 
goal,  whereas  the  present  method  uses  both  an  NMC  goal  and 
an  SDO  goal,  and  that  WARS  can  accomodate  non-cannable  and 


Indentured  items. 


The  basic  structure  of  the  VARS  mode]  is  depicted  in 
Figure  1-1.  For  simplicity  the  discussion  will  he  limited 
to  aircraft  spare  parts,  although  WARS  is  applicable  to  all 
types  of  WRM.  The  required  innut  data  in  block  1  includes: 
General 

Period  covered  (days) 

No.  of  assigned  aircraft 

Dally  Flying  Hour  Program  (DFHP) 

Maximum  expected  NMC  goal 

Individual  Item  Data 
National  Stock  Number  (NSN) 

Parent  item  NSN 

Hierachy  position  on  aircraft 

Type  I-an  SP.U  assigned  to  an  LRU 
Type  Il-an  LRU  assigned  to  a  module 
Type  Ill-an  LRU  or  module  assigned  to  an  engine 
Type  IV-an  LRU  or  engine  assigned  to  an  aircraft 
Kind  of  item 

LRU,  SRU,  engine,  module 

reparable, irreparable 

cannable,non-cannable 

quantity  per  higher  assembly  (OPHA) 

quantity  per  aircraft  (QPA) 

unit  cost 

Base  Repair  Cycle  (BP.C) 

Depot  Repair  Cycle  (DRC) 
transportation  time  to  depot 
total  demand  per  hundred  flying  hours 
Base  Repair  Rate  (BRR) 

Base  Condemnation  Rate  (BCR) 

Depot  Condemnation  Rate  (DCR)  (Ref  9:  8) 


block  2  to  calculate  the  max  in  tin  expected  number  of  each 
item  "that  have  been  removed  from  the  aircraft  due  to 
failures  and  have  not  yet  been  returned  to  a  serviceable 
condition  available  for  replacement  onto  the  aircraft." 
This  quantity  is  called  the  "pipeline  quantity"  or  "MU"  and 
is  the  expected  value  on  which  the  probability  distribution 
of  demands  is  based.  (Ref  9:  10)  For  example,  if  the  sup¬ 
port  period  under  consideration  is  ten  days  and  the  kit  con¬ 
tains  two  parts,  the  PIPE  subroutine  night  calculate  the 
foil owi ng  data: 

EXPECTED 

DAY 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

Table  1-2 


NO.  OF  PARTS 
PART  A 
0 

.3 
.7 
1  . 1 
.8 
.7 
.7 
.6 
.6 
.  5 


UNAVAILABLE 
PART  R 
.1 
.4 
1.3 
2.7 
3  .  1 
2.6 
2.2 
2.2 
2 . 1 
2.1 


Note:  the  decrease  around  day  four  could  again  be  caused  by 

the  arrival  of  repair  capability  or  a  decrease  in  the  PFHP. 


The  naxlmis  expected  number  of  Part  A  unavailable  is 
1.1  which  occurs  on  day  four,  while  the  maximum  exnected 
numberr  of  Part  P»  unavailable  is  3.1  which  occurs  on  day 
five.  So  the  pipeline  quantities  or  MUs  for  Parts  A  and  B 
are  1.1  and  3.1  respectively. 

Block  3  uses  the  MU  values  calculated  in  PIPE  to  deter¬ 
mine  the  expected  number  of  KMC  aircraft  using  the  "ASSESS" 
assessment  subroutine,  based  on  the  present  stock  in  the 
kit.  (Ref  :  11)  The  PIPE  and  ASSESS  subroutines  will  be 
explained  in  greater  detail  later. 

Block  A  checks  to  see  if  the  E ( N M C ]  goal  set  by  the 
operator  is  met.  If  so,  the  kit  stock  levels  and  the  E[NMC) 
figure  are  output,  otherwise  the  MA  process  is  performed  to 
determine  which  item  should  be  added  to  the  kit  to  obtain 
the  greatest  decrease  in  EfNMC]  per  dollar  spent  (block  5). 
Once  this  is  accomplished,  the  kit  is  re-evaluated  by  ASSESS 
(block  3)  and  the  process  is  repeated. 

PIPE  calculates  the  daily  expected  number  of  each  item 
removed  from  aircraft  due  to  failure  and  not  yet  service¬ 
able.  MU  is  then  set  equal  to  the  maximum  of  the  daily 
quantities  for  the  desired  period  (Ref  9).  Figure  1-2  depi- 
cits  the  flow  of  failed  parts  through  the  "pipeline.” 
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Parts  Pipeline 

Figure  1-2. 


Suppose  a  certain  part  fails  on  an  aircraft.  It  is 
removed  and  enters  the  pipeline.  Usually  the  part  will  be 
repaired  at  base  level,  however  a  certain  percentage  will 
require  depot  level  maintenance.  The  percentages  that  go  to 
the  base  and  to  the  depot  are  called  the  Base  Repair  Rate 
(BRR)  and  the  Depot  Repair  Fate  (PRR)  respectively.  At  both 
the  base  and  the  depot  a  certain  percentage  will  be  con¬ 
demned.  These  percentages  are  known  as  the  Base  Condemna¬ 
tion  Rate  (BCR)  and  the  Depot  Condemnation  Rate  (DCR) 
respectively.  The  time  required  to  repair  a  particular  item 
at  base  level  is  called  the  Base  Repair  Cycle  (BRC)  and 
similarly  a  Depot  Repair  Cycle  (DRC)  at  the  depot  level. 
Based  on  the  DFHP  and  the  above  input  parameters,  PIPE  cal¬ 


culates  MU  for  each  i  • em 


The  cornerstone  of  the  VA  process  of  optinization  is 


the  ability  to  calculate  the  E  f  N  M  C 1  for  any  stock  level  in  a 
kit.  The  VARS  model  performs  this  task  using  the  subrou¬ 
tine  ASSESS.  Once  the  maximum  pipeline  quantities  are 
determined  for  each  part,  this  information  is  used  by  the 
ASSESS  subroutine  to  calculate  the  maximum  E[NMC)  for  a 
given  stock  level  using  the  following  expected  value  for¬ 
mula: 

EfNMC]*0xP[0  NMC  a/c]  +  lxP[l  KMC  a/c]  +  2xP[2  KMC  a/c]  + 

3xP  [  3  NMC  a/c]  +  ...  +  UExPfUE  NMC  a/c]  (1) 

where  UE  is  the  number  of  aircraft  assigned  to  the  unit  (Ref 
9:  16).  Computing  the  above  probabilities  is  extremely  com¬ 
plex  because  of  the  possibility  of  some  parts  being  canni¬ 
balized  from  other  aircraft  and  the  indenture  relationship 
between  different  parts. 

Figure  1-3  shows  some  indenture  relationships  in  a  sim¬ 
ple  component  hierarchy  for  a  fictitious  aircraft.  Note 
that  SP.U  5  and  SRU  6  are  component  parts  of  LRU  S  which  is 
itself  a  component  part  of  module  1.  Also  note  that  the 
type  of  item  is  given  on  the  line  above  the  item.  For  exam¬ 
ple,  SRU  5  and  SRU  6  are  type  I  items  while  LRU  5  is  a  type 
II  item.  Recall  that  type  designations  merely  reflect  the 
indenture  level  of  the  part. 


When  discussing  indenture  relationships  it  is  important 


to  understand  that  if  an  SRI’  such  as  SRU  3  has  failed  and  no 
replacement  is  available,  its  parent  parts  up  the  tree  (i.Rt: 


3  and  engine  1)  also  fail.  The  converse  is  not  true.  If  a 
parent  part  such  as  engine  1  fails,  it  has  not  necessarily 
failed  because  of  a  failure  of  one  of  its  component  parts. 
There  are  numerous  pieces  in  an  engine  assembly  itself 
besides  the  lover  indentured  items  which  could  fail.  Said 
another  way,  all  the  parts  that  make  up  a  parent  item  are 

not  identified  as  separate  components  and  shown  on  a  parts 
hierarchy. 
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P.efore  discussing  the  computational  methodology  some 
additional  terms  need  to  be  defined  (Ref  10): 

higher  assembly  (HA)-refers  to  an  item  that  contains 
parts  that  can  only  be  reached  by  taking  this  item 
apart,  sometimes  called  a  parent  with  respect  to 
immediately  lower  parts. 

quantity  per  higher  assembly  (OPHA)-the  quantity  of  a 
given  part  contained  within  each  of  its  parent 
parts. 

quantity  per  aircraft  (QPA)-the  quantity  of  a  given 
part  contained  on  one  aircraft. 

backorder-an  unfilled  demand  or  shortage  at  base  level. 

perfectly  cannibal izable  part  (PCP)-a  part  that  can  be 
cannibalized  from  one  I! A  to  another  and  which  has 
the  same  OPHA  on  all  HAs  containing  the  item. 

non-cannibal izable  part  (NCP)-  a  part  that  can  not  be 
removed  in  a  serviceable  condition  from  one  air¬ 
craft  or  11A  to  prevent  a  backorder  on  another  air¬ 
craft.  While  it  may  be  physically  possible  to 
cannibalize  such  an  item,  the  item  may  be  desig¬ 
nated  as  a  NCP  because  the  item  would  be  damaged 
when  it  is  removed  or  the  time  to  cannibalize  it 
is  considered  excessive. 

imperfectly  c anni ba 1 i za bl e  part  (MCP)-a  part  that  can 
be  cannibalized  from  one  HA  to  another  and  which 
does  not  have  the  same  OPHA  on  all  HAs  containing 
the  item.  Also  called  myopically  cannibal i za bl e 
par t . 

national  stock  number  (NSN)-a  unique  number  by  which 
each  part  procured  by  the  federal  government  is 
identified . 

unique  item-an  NSN  that  appears  in  only  one  position  in 
the  aircraft  parts  hierarchy. 

non-unique  item-an  NSN  that  is  used  in  more  than  one 
hierarchial  position  on  an  aircraft. 


A  comprehensive  list  of  terms  is  given  in  the  glossary  in 


appendix  A 


In  describing  the  assessment  procedure  the  steps  are 
presented  in  the  reverse  of  the  order  of  occurrence.  This 
is  done  to  continuously  provide  the  reader  with  the 
rationale  for  each  successive  step.  A  detailed  account  of 
the  variables  and  formulas  used  in  the  assessment  routine  is 
included  in  appendix  R.  This  general  description  will  start 
with  the  final  answer  and  work  backward  through  lower  and 
lower  indenture  levels  as  additional  information  is  needed. 
In  the  actual  routine  the  lowest  level  is  evaluated  first 
and  the  information  obtained  from  it  used  in  each  successive 
level  up  the  tree  until  the  required  Information  is  avail¬ 
able  at  the  top  level  to  compute  E[NMC).  The  procedure  will 
be  explained  for  the  top  indenture  level  but  the  other  lev¬ 
els  are  evaluated  in  a  similar  manner  (Ref  9,11). 

E(NMC]  for  a  given  stock  level  is  found  using  the  fol¬ 
lowing  steps: 

1.  If  U£  is  the  number  of  aircraft  in  the  unit  under  con¬ 
sideration,  then  using  the  expected  value  formula  for 
discrete  random  variables  (Ref  12:  9f>): 

E  f  NMC  ]  “OxP  ( 0  NMC  a/c]  +  lxP(l  NMC  a/c]  +  2xP[2  NMC  a/c)  + 

3 x P [ 3  NMC  a/c]  +  ...  +  UExPfUE  NMC  a/c]  (1) 

2.  The  probabilities  that  a  certain  number  of  aircraft  are 
NMC  are  merely  the  probability  function  (PF)  of  the  number 


of  aircraft  NMC  and  can  easily  be  calculated  if  the  cun  illa¬ 
tive  distribution  function  (CDF)  is  known  (Ref  13:  f>0). 
Therefore  if  the  CDF  is  known  the  TF  can  be  calculated  and 
the  individual  probabilities  used  in  the  right  side  of  eciua- 
t  i  on  1  . 

3.  The  CDF  needed  in  step  2  is  obtained  by  computing: 

p ( A  a/c  NMC<x)  =  P (#  a/c  KMC  for  MCP  or  MCP<x)  x 

P (»  a/c  NMC  for  PCP<x)  (2) 

For  all  x*0 , 1 , 2 ,  .  .  .  , UE 

The  number  of  aircraft  NMC  is  the  maximum  of  the  number  of 
aircraft  NMC  for  PCP  and  the  number  of  aircraft  NMC  for  NCP 
or  MCP.  So  symbolically: 

P ( Ma x [ Y , Z ] <c )  -  P ( Y< c  and  Z<c)  -  P(Y<c)  x  P(Z<c) 

In  effect  this  is  multiplying  two  independent  CDFs  by  ele¬ 
ments  to  arrive  at  a  resultant  CDF.  The  CDFs  are  combined 
in  this  manner  because  the  distribution  of  the  maximum  of 
the  two  random  variables  (max[Y,Z])  is  required  and  not  the 
sum  (Y+Z)  (Ref  12:  253,  14:  310).  This  is  true  because  can¬ 
nibalization  can  be  used  to  consolidate  the  failures  onto 
fewer  airframes.  Consider  an  example  combining  the  follow¬ 


ing  two  CDFs : 


V*  ■- 


CDF  1 


0  1 


P ( X<x )  .1  .5  1.0 


CDF  2 


0  1 


P(X<x)  .3  .7  1.0 


Multiplying  elements  yields  the  following: 


RESULTANT  CDF 


P ( X<x )  .03  .35  1.0 


Converting  this  CDF  to  a  PF : 


RESULTANT  PF 


P(X-x)  .03  .32  .65 


The  sane  answer  could  be  obtained  by  using  the  following 
PFs  , 
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X 


P(X«x)  .1  .4  .5 

PF  2 

x  0  12 

P  (  X  =  x  )  .3  .4  .3 

and  performing  the  following  calculations: 

P(X*0)  ■  .1  x  .3  *  .03 

P(X**1)  *  .1  x  .4  +  .4  x  .3  +  .4  x  .4  =  .32 

P(X*=2)  *  .  1  x  .  3  +  .  4  x  .  3  +  .  5  x  .  3 

+  .5  x  .3  +  .5  x  .4  *  .  f>  5 

These  results  agree  with  the  resultant  PF  above  as  expected. 

The  CDFs  are  independent  rather  than  mutually  exclusive 
because  an  aircraft  could  he  broken  for  several  different 
types  of  parts  at  the  same  tine.  For  example,  an  aircraft 
could  initially  be  broken  for  an  NCP,  and  subseouently  have 
a  PCP  cannibalized  from  it.  Thus  the  aircraft  would  be  bro¬ 
ken  for  two  different  types  of  parts. 

4.  At  this  point  the  procedure  splits  into  two  branches. 
We  will  first  trace  through  the  procedure  for  determining 
the  last  term  in  equation  2  and  return  later  to  explain  how 


the  other  term  on  the  right  hand  side  of  the  equation  is 


obtained.  The  stens  in  the  first  branch  will  be  numbered 


'5a",  “6a",  etc.  while  the  second  branch  will  be  numbered 


'  5  b " ,  “6b",  etc. 


5a.  Tracing  the  first  branch: 


P  ( a/c  NM  C  for  PCP^x)  *■  P  ( 1‘  a/c  KMC  for  PCP  l^x)  x 


P(0  a/c  KMC  for  TCP  2^x)  x 


P(*  a/c  KMC  for  PCP  last/x) 


For  all  x*0 , 1  ,  2  ,  . . . , UE 


The  CDFs  on  the  right  side  of  the  equation  are  again 


independent,  and  it  is  the  maximum  of  the  random  variables 


which  is  of  interest  using  the  sane  logic  as  in  step  3 


above  . 


At  this  point  only  the  top  indenture  level  narts  are  con¬ 


sidered  and  only  PCPs  . 


6a.  Each  of  the  probability  expressions  in  equation  3  can 


be  calculated  by  the  formula 


P(«  a/c  KMC  for  PCP  i  ^x) 


x 

2 


y  -u 


Again  for  all  x*0,l,2,...UE 


If  it  is  assume^  that  the  number  of  NMC  aircraft  for  any  I’ CP 


is  Poisson  distributed  and  the  rtoan  is  known  for  each  part. 
Since  the  CPF  in  equation  four  is  of  finite  leng  th,  WARS 
adds  the  probabilities  in  the  far  tail  of  the  Poisson  dis¬ 
tribution  into  the  last  interval  used.  ASSESS  currently 
assumes  a  Poisson  distribution  but  it  will  be  able  to  use 
the  Negative  P.inomial  distribution  in  a  later  version. 

7a.  The  "u"  or  TOTMU  values  for  each  PCP  used  in  eouation  4 
are  determined  by  using  the  formula: 

TOTMU  “  AMU  +  AW P  (5) 

where  AMU  is  the  "actual  MU"  or  expected  pineline  quantity 
for  the  part  calculated  by  the  PIPE  subroutine  and  AWP  is 
the  number  "awaiting  parts"  or  unavailable  due  to  a  lack  of 
lower  indentured  parts.  The  AWP  figures  for  each  nart  are 
the  only  results  of  the  calculations  for  a  lower  indenture 
level  used  in  the  next  higher  level.  All  items  with  no 
subassemblies  have  an  AWP  figure  of  zero.  For  the  top  level 
there  is  only  one  AWP  figure  and  is  called  E [ N  M  C  ]  ,  the 
number  of  aircraft  awaiting  lower  indentured  parts. 

5b.  Returning  to  step  4  and  tracing  the  branch  for  the 
other  term  on  the  right  hand  side  of  equation  2: 
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P(0  a/c  KMC  for  NCP  or  MCP^x) 


N  M  C  for  NCP  or 


P  (KMC 


for  NCP  or 


UE-n 

f'crl]  (6) 


A  CDF  is  obtained  by  using  this  binomial  formula  for 
x“0,l,2,...,UE.  For  the  top  indenture  level  there  is  only 
one  HA,  the  aircraft  itself,  and  hence  one  CDF.  For  lower 
levels  there  will  be  a  CDF  for  each  HA  in  that  level  and  UE 
in  the  equation  will  be  replaced  by  the  quantity  of  each  HA 
in  the  fleet. 


6b. 


P ( NM C  for  NCP  or  MCP)  «  l-[p(all  NCP  avail.)  x 

P  (  al  1  MCP  avail  .7] 


(7) 


7b. 


P(all  NCP  avail.)  =  P(lst  NCP  avail.)  x 

P(2nd  NCP  avail.)  x 


P ( 1  a  s  t  NCP  avail.) 


(8) 


2  U 


■fcp.W  .Vv'A*'/V  V  V  v  V  V  V  v  v  V  v 


1  -S  A  -Y  A  . 


P ( a  1 1  MCP  avail.)  =  P(lst  MCP  avail.)  x 


P  (  2  n  d  MCP  avail.)  x 


P(last  MCP  avail.) 


(9) 


8b.  In  equation  8  if  all  parts  have  unique  QPHAs: 


P( ith  NCP  avail . 


expected  shortages(i) 


1- 


quantity  of  part  i  in  fleet 


op: 

T_ 


1 A  (  i  ) 
(10) 


with 

expected 


8 1  ock  1 evel  +  qty  on  fleet 
shortaRes*^^  (x-stock  level)  P(x  demands) 
x*  s  to  ck  1  eve  1 


(11) 


and  P(x  demands)  determined  from  the  total  pipeline  quantity 
in  equation  5.  The  bracketed  expression  in  equation  ten 
will  never  be  negative  because  the  expected  shortage  can 
never  exceed  the  quantity  in  the  fleet. 

If  there  are  some  multiple  QPHA  parts,  equation  10  becomes 
slightly  more  complex. 

In  equation  9: 


expected  #  of  HAs  missing  items 

P  (  i  t  h  'TCP  avail.)-  1-  _ _ _  (12) 

total  If  of  11  As  o"f  part  I 


The  numerator  in  equation  12  can  he  found  if  the  PF  is 
known.  The  PF  can  be  calculated  since  the  PF  of  the  number 
of  demands  is  known  and  the  number  and  PPHA  of  each  HA  is 
known.  ASSESS  assumes  all  MCPs  are  cannibalized  first  from 
the  HAs  with  the  greatest  quantity  inside. 

The  PFs  of  the  number  of  demands  used  in  equations  11 
and  12  all  need  the  TOTH Us  from  equation  5.  Therefore  the 
AUP  result,  or  E[NMC),  of  the  top  level  depends  solely  on 
the  AMUs  provided  by  the  PIPE  subroutine  and  the  AWT 
results  of  the  next  lower  level. 

The  preceding  assessment  methodology  makes  the  follow¬ 
ing  approximations  (Ref  15,16): 

1.  The  probability  distributions  are  reset  to  Poisson 
from  one  indenture  level  to  another  as  explained 
be  1 ow . 

2.  The  backorders  of  a  nonunique  item  have  a  Poisson 
distribution  in  each  hierarchy  position  with  that 
stock  number. 

3.  The  probability  distribution  of  HAs  awaiting  parts 
due  to  shortages  of  N'CP  or  MCP  is  binomial. 

It  is  the  first  assumption  which  is  of  interest.  For 
each  part  i,  the  WARS  assessment  routine  calculates  the  PF 
(actually  the  CDF)  of  the  number  of  parts  AWP  and  uses  it  to 
compute  a  mean  AWP  (frequently  called  Expected  Back  Order 


(EPt')).  Adding  the  mean  of  this  A  UP  distribution  to  thr 
nean  of  the  pipeline  ouantitv  for  part  i  in  equation  5 
results  in  a  combined  mean  which  is  used  to  generate  a  Pois¬ 
son  distribution  for  its  parent  part.  This  Poisson  distri¬ 
bution  is  used  as  an  approximation  of  the  true  distribution. 
The  demands  of  the  subassemblies  of  part  i  are  Poisson  but 
when  the  stock  level  of  the  subassenblifes  are  taken  into 
account  the  backorder  distribution  may  not  be  Poisson.  Com¬ 
bining  subassembly  backorder  distributions  results  in  an  AWP 
distribution  for  part  i  that  may  not  be  Poisson. 

The  correct  way  of  combining  the  A'.IP  and  the  AMU  PFs  to 
get  a  TOTMU  PF  is  as  follows: 

P(TOTMU-O)  -  P ( AWP*0 )xP(AMU*0) 

P  ( TOTMU  *1  )  -  P(  Ab’P*0)xP(  AMU-1  )  +  P  (  AVP-1  )xP  ( AMC-0  ) 

P(T0TMl>2)  -  P(AWP=0)xP(AMU  =  2  )  +  P ( AVP=1 ) x P ( AMU-1  ) 

+  P(AWP-2)xP(AMU-0)  (13) 


This  Po is sonSappr oximat ion  of  the  number  of  demands  of 
a  parent  part,  when  carried  through  several  indenture  lev¬ 
els,  may  induce  a  significant  error  in  the  final  E [ N M C 1  fig¬ 
ure  and  could  therefore  degrade  the  ability  of  ASSESS  to 
evaluate  the  performance  of  a  kit.  Since  the  optimal  kit 


building  capability  of  WARS  hinges  on  tbe  ability  to  accu- 
ritelv  evaluate  kit  performance,  the  effect  of  this  approxi¬ 
mation  is  obviously  an  important  area  to  investigate. 

AFLC/XRS  has  made  a  cursory  check  of  most  of  the 
approximations  used  in  the  assessment  routine  hut  the  check 
was  only  made  for  one  relatively  small  parts  hierarchy  (Ref 
15).  The  Poisson  approximation  was  always  evaluated  in  con¬ 
junction  with  other  approximations  and  was  never  evaluated 
separately . 

The  objective  of  this  study  is  not  to  evaluate  how  well 
WARS  builds  optimal  WP.M  packages  nor  is  it  to  evaluate  how 
well  its  assessment  routine  can  determine  a  kit's  perfor¬ 
mance.  Rather,  the  purpose  of  this  thesis  is  to  investigate 
the  effect  of  resetting  the  demand  probability  distribution 
of  every  II A  to  a  Poisson  distribution  at  each  indenture 
level.  The  approximation's  effect  on  the  WARS  assessment 
results  will  be  evaluated  for  a  variety  of  conditions  to 
determine  its  sensitive  parameters. 

Scope 

Although  several  approximations  are  made  in  the  assess¬ 
ment  routine,  the  Poisson  approximation  will  be  the  only  one 
investigated.  In  order  to  eliminate  the  effects  of  the 
other  approximations  on  the  results,  this  study  will  he  lim¬ 
ited  to  fully  cannibalizable,  unioue  parts  with  unique 


QPHAs.  If  non-cannibalizahle  parts  were  allowed,  the  effect 
of  approximation  3  would  also  be  included  in  the  results. 
If  non- unique  parts  were  allowed,  approximation  2  would 
affect  the  results.  Likewise,  if  narts  with  several  dif¬ 
ferent  Q  Pit  As  were  allowed,  the  results  would  apain  be 
affected  by  approximation  3. 

The  WARS  model  builds  optimal  packages  for  WRSK,  BLSS, 
and  OWRfl.  Since  differences  between  these  kits  are  incor¬ 
porated  into  the  part  pipeline  quantities,  our  study  will 
encompass  all  types  of  WRM. 

Approach  and  Presentat ion 

To  evaluate  the  effect  of  the  approximation  two  com¬ 
puter  models  were  developed.  One  model,  called  APPROX, 
evaluated  a  package  of  parts  using  the  WARS  procedure  with 
the  Poisson  approximation.  The  other  model,  called  EXAOT, 
evaluated  a  package  of  parts  using  exact  probability  calcu¬ 
lations.  Both  models  output  the  demand  probability  distri¬ 
butions  at  each  hierarchial  position  as  well  as  the  final 
EfNNC]  figure. 

To  determine  the  effect  of  the  approximation,  the  same 
data  was  input  to  both  models  and  the  differences  noted.  By 
constructing  different  aircraft  hierarchies  and  using  dif¬ 
ferent  part  data,  sensitivity  to  the  following  characteris¬ 


tics  was  evaluated: 


depth  of  indenture  of  the  ports  hierarchy 


-  S  A  to  ! ! A  ratio  for  higher  levels 

-  stock  level  of  each  part 

-  QP1IA  of  each  part 

All  numerical  results  are  included  in  appendices  E 
through  G,  and  summaries  of  the  results  are  presented  in 
chapter  three  in  tabular  form. 
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In  order  to  quickly  and  accurate! v  calculate  the  effect 
of  the  Poisson  approxination  for  different  pipeline  quanti¬ 
ties,  stock  levels,  and  hierarchlal  structures,  two  computer 
proprans  were  developed.  The  first  program,  called  APPROX, 
uses  the  same  computations  and  Poisson  approximation  as  the 
WARS  ASSESS  subroutine.  The  other  program,  called  EXACT, 
uses  similar  computations  but  calculates  the  exact  probabil¬ 
ities  by  not  using  the  approximation.  Both  programs  evalu¬ 
ate  backorder  PFs  for  each  part  in  the  aircraft  parts 
hierarchy  assuming  all  parts  are  c ann ibal i za b le ,  all  parts 
are  unique,  and  no  parts  have  multiple  OPHAs.  The  two  pro¬ 
grams  are  described  below. 


APPROX  Computer  Program 


APPROX  evaluates  a  certain  parts  hierarchy  one  part  at 
a  time  beginning  with  the  lowest  indentured  parts,  evaluat¬ 
ing  an  entire  indenture  level  before  moving  up  to  the  next 
level.  It  does  this  by  storing  all  data  and  probability 
distributions  in  a  three-dimensional  array  called  "A" .  Fig¬ 
ure  2-1  illustrates  the  array  format  used  to  store  each 
part's  information.  This  figure  will  be  used  later  to  help 
explain  the  logical  steps  of  APPROX.  The  actual  APPROX  pro¬ 
gram  and  a  listing  of  computer  variables  used  in  the  program 


are  included  in  appendix  C 


The  following  is  a  list  of  definitions  of  the  acronvns 


used  in  figure  2-1: 


A  -  t hr ee-d ine ns io na 1  data  and  probability  array 


AMU  -  A  variable  denoting  the  basic  mean  pipeline  quantity 
for  an  item,  excluding  AWP.  Also  called  "actual  mu"  or 


APPROX  -  program  which  calculates  approximate  F.  [  N  M  C  ]  • 


AUP  -  A  variable  denoting  the  expected  number  of  items 
“awaiting  narts." 


BO  -  backorder. 


CDF  -  Cumulative  Distribution  Function 


DISTL  -  desired  length  of  probability  distributions  in  array 


EBO  -  expected  number  of  backorder. 


EXACT  -  program  which  calculates  exact  EfNMC] 


HA  -  Higher  Assembly. 


MU  -  pipeline  quantity  from  which  the  Poisson  distribution 


is  generated.  Used  i n t e r c h a n ge ah  1 v  with  AMU. 


PF  -  Probability  Function 


PUP  -  position  of  a  suhassonhlv  under  a  parent  iron,  for 
example,  if  a  parent  item  has  two  stihassenhlies(SA), 
the  first  S  A  will  have  a  PUP  of  1  and  the  second  5 A 
will  have  a  TUP  of  2. 

QPHA  -  Quantity  Per  Higher  Assembly.  This  is  the  number  of 
each  item  on  its  respective  higher  assembly. 

SA  -  Subassembly. 

SL  -  stock  level  for  a  specific  part. 

TOTMU  -  A  variable  denoting  the  total  expected  pipeline 
quantity  for  an  item.  It  is  equal  to  AMU  +  AWP. 

Type  (of  an  item)  ~  The  type  of  an  item  refers  to  its  inden¬ 
ture  relationship  as  follows: 

4  **  An  LRU  or  engine  attached  directly  to  the  aircraft. 

3  *  An  LRU  or  module  attached  to  a  type  4  item. 

2  *  An  LRU  attached  to  a  module. 

1  *  An  SRU  attached  to  any  LRU. 


all  suhassenblies  of  this  part  into  an  A'.'P  CPF  for  this 
pa  r  t  . 

2.  The  Ah' P  CDF  (line  7)  is  converted  to  an  AVP  PF 
(line  S )  . 

3.  The  expected  nunher  of  parts  AWP  is  calculated  fron 
the  AWP  PF  (line  8)  and  stored  in  line  9  column  1. 

4.  Th  expected  number  of  narts  A'-’P  from  step  3  is 
added  to  the  expected  pineline  quantity  for  this  part  (AMU) 
from  line  1,  column  4,  and  the  result  (TOTMl')  is  stored  in 
line  9,  column  2. 

5.  TOTMU  from  step  4  i6  used  as  the  parameter  to  gen¬ 
erate  a  Poisson  TOTMU  PF  which  is  stored  in  line  10, 

6.  The  TOTMU  PF  is  converted  to  a  TOTM]  CDF  (line  11). 

7.  The  TOTMU  CDF  from  line  11  is  adjusted  for  stock 
level  and  QPHA  to  yield  a  BO  CDF  for  this  part  (line  12). 

8.  The  BO  CDF  (line  12)  is  converted  to  a  BO  PF  (line 

13)  . 

9.  The  variance,  mean,  and  variance  to  mean  ratio  for 
the  BO  PF  are  calculated  and  stored  in  line  14,  columns  0 
through  2  respectively. 

10.  The  BO  CDF  (line  12)  is  transferred  to  its  parent 
part  to  be  used  as  a  subassembly  BO  CDF  in  step  1. 


o 


n 

These  sane  10  steps  arc  (’one  for  each  part,  but  thev  must 
be  done  for  the  tvpe  I  (most  deeply  indentured)  parts  first. 
After  all  type  I  parts  have  been  processed,  the  type  II 
parts  are  considered,  and  so  on  through  tbe  higher  part 
types.  The  highest  part  type  in  WARS  is  type  IV,  but  for 
the  purposes  of  the  APPROX  program  the  aircraft  itself  is 
considered  a  type  V  part  and  is  processed  last. 

To  better  understand  the  10  steps  above,  an  example  is 
presented.  Consider  the  parts  hierarchy  shown  below  and  its 
associated  parts  data: 


PARTS 

DATA 

Part 

OPHA 

MU 

Stock 

H 

1 

1 

2 

A 

1 

1 

0 

R 

2 

2 

1 

C 

2 

1 

0 

C  . 


The  computations  start  with  the  >. ,  I  parts  A,  B,  and 
For  part  A  there  are  no  subassemblies  so  the  subassembly 


2 


& 


-■ .  w:  v. 


vy.'.'.v. 


x 


0 


3  4 


P(^FO^x)  1.000  1.000  1.000  1.000  1.000  .  .  . 


Multiplying  then  colunruise 

In  step  1 , 

yields  an  AUP  CDF  of 

X 

0 

1 

2 

3 

4 

P ( # AWP^x ) 

1 .000 

1 .000 

1 .000 

1 .000 

1.000  .  .  . 

In  step  2  the 

AUP  CDF 

Is  converted  to 

the  following  AUP  PF: 

X 

0 

1 

2 

3 

4 

P(/»AWP^x) 

1 .000 

0.000 

0.000 

0.000 

0.000  .  .  . 

Step  3,  the  expected 

number 

of  parts  AVP  is 

calculated  to  be 

0. 

Step  4,  TOTMU 

-  AMU  + 

AUP  - 

1+0-1. 

Step  5,  using 

TOTMU  » 

1  ,  the 

f ol lowing 

Poisson  TOTMU  PF  is 

generated : 

X 

0 

1 

2 

3 

4 

P(//de:nands*x) 

.3679 

.  3679 

.1839 

.0613 

.0153  .  .  . 

Sten  6,  the  above  PF 

Is  converted  to  a 

TOTMU 

CDF  : 

X 

0 

1 

2 

3 

4 

P(#d emands^x ) 

.  3679 

.  7  358 

.9197 

.9810 

.9963  .  .  . 

Step  7 ,  since 

stock  * 

0  and 

QPH A  =  1, 

the  F0 

CDF  is  identi- 

cal  to  the  TOTMU  CDF: 

X 

0 

1 

2 

3 

4 

p( ftno^x) 

.3679 

.  7  358 

.9197 

.9810 

.9963  .  .  . 

3  7 


Step  8  ,  the  BO  TF  is  then: 


X 

0  1 

2  3 

4 

P(/'RO«x) 

.3679  .3570 

.1839  ,0613 

.0153  .  .  . 

Step  9 

Variance  “1 

Mean  or  EBO  *  1 

Variance  to  mean 

ratio  *  1 

Sten  10,  the 

BO  CDF  calculated 

in  step  7  will 

he  the  distri- 

btition  used 

in  the  calculations 

for  part  H . 

Steps  eight  and 

nine  only  need  to  be  accomplished  for  the  aircraft  itself, 
but  are  performed  for  each  part  for  analysis  reasons. 


For  part  B,  the  results  of  the  steps  are: 
Step  1 


X 

0 

1 

2 

3 

4 

P(/'AWP^x) 

1.000 

1.000 

1 .000 

1 .000 

1.000  .  .  . 

Step  2 

X 

0 

1 

O 

4. 

3 

4 

P(#AVP*x) 

1 .000 

0 .000 

0.000 

0.000 

0.000  .  .  . 

Step  3,  AW P  » 

Step  4,  TOTMU 

0  . 

-  AMO  + 

AW  P  =■  2 

+  0  *  2 

• 

Step  5 

X 

0 

1 

2 

3 

4 

P(#demands*x) 

.1353 

.  2707 

.  2  702 

.1805 

.0902  .  .  . 

Step  6 


X 

0 

1 

2 

3 

4 

P ( ^d  ema  nd  s^x  ) 

.1353 

4  06  6 

.  6767 

.  8  57  1 

.  947  3 

Step  7,  since 

the  stock 

level 

-land 

the  QPH A 

-  2, 

CDF  is 

X 

0 

1 

2 

3 

4 

P(#B0^x) 

.4060 

857  1 

.  9  R  3  4 

.  9989 

.9999 

Step  8 

X 

0 

1 

2 

3 

4 

P(#B0-x) 

.4060 

.4511 

.1263 

.0155 

.0011 

Step  9  Variance  *  .5442 

Mean  or  EBO  -  .7546 

Variance  to  mean  ratio  *  .7211 

Sten  10,  the  F0  CDF  from  steo  7  will  be  used  in  calculations 
for  part  11. 

For  Part  C 

Steps  1  through  3  again  yield  the  same  results  as  above. 

Step  4,  TOTMU  -  AMU  +  AWP  *1+0-1. 

Step  5 

x  0  12  3  4 

P(Nenands-x)  .3679  .  367  9  .  1839  .  061  3  .01  53  .  .  . 


P  (  #  dema  nd  s^!x  ) 

.  3  670 

.  7  3  5  fi 

.9197 

.9810 

.9963  .  .  . 

Step  7,  stock 

level  ■ 

0  and 

0  PM  A  -  2, 

s  0 

X 

0 

1 

2 

3 

4 

Pf/'RO^x) 

.  3679 

.9197 

.  9963 

.  9999 

1.000  .  .  . 

Step  8 

X 

0 

1 

2 

3 

4 

P(  /'RO-x) 

.  3679 

.5518 

.  0766 

.0036 

.0001  .  .  . 

Step  9  Variance  *  .3790 

Mean  or  EBO  *  .7162 
Variance  to  mean  ratio  =  .5292 
Step  10,  again  the  BO  CDF  from  step  7  will  be  used  in  tl 
calculations  for  part  H. 

Now,  type  II  parts  are  considered. 

For  part  H 

Step  1,  using  the  subassembly  P0  CDFs  previously  calculated 


X 

0 

1 

o 

A. 

3 

4 

P  (  #  A 

BO^X) 

.  3679 

.  7  3  58 

.9199 

.9810 

.9963  . 

P(#B 

BO^x) 

.4060 

.  8571 

.9834 

.9989 

.9999  . 

P(i»C 

BO^x) 

.  3679 

.  9897 

.  9963 

.9999 

1.000  . 

P(0H 

AWP^x) 

.0549 

.  5900 

.9012 

.  9799 

.9963  . 

S  t  on  2 


X 

0 

1 

4. 

3 

P(/fI!  AVJP=x) 

Step  3,  A VP  « 

.0  549 

1.4885. 

.  52  5  1 

.3212 

.0787 

.0 

Step  4,  TOTMU 

■=  AMlT  + 

AWP  »  1 

.0  +  1  . 

4885  =  2. 

4885 

Step  5 

X 

0 

1 

2 

3 

P(^denands*x) 

.0830 

.  2066 

.  2  571 

.2133 

.  1 

Step  6 

X 

0 

1 

2 

3 

P( flderoands^x) 

.0830 

.  2897 

.  5468 

.7600 

.89 

Step  7 

X 

0 

1 

3 

P(/>BO^x) 

.  5468 

.2133 

.1327 

.0660 

.02 

Step  8 

x 

0 

1 

2 

3 

P(0P>O*x) 

.  5468 

.2133 

.1327 

.0660 

.0 

Step  9  Variance  =  1.4466 

Mean  or  EBO  *  .8612 

Variance  to  mean  ratio  =  1.6799 

In  this  example  there  are  no  type  III,  IV,  o 
so  the  EBO  of  part  H  is  the  desired  result. 


EXACT  Conn u ter  Pronran 

The  EXACT  computer  procram  is  also  included  in  appendix 
C.  This  procram  is  similar  to  the  APPROX  program  and  uses 
the  same  computer  variables  previously  listed  for  APPROX. 

Figure  2-2  shows  the  array  format  for  EXACT.  Note  that 
it  differs  from  the  APPROX  format  only  on  line  9,  where  the 
AWP  and  TOTMU  figures  are  replaced  by  the  Ax,.l’  PF . 


Format  for  Each  Fart  in  Array  A 


ERO 

Va  r 

Mean 

Figure  2-2 

The  computational  steps  used  in  F.XACT  are  similar  to 
those  in  APPROX  except  that  steps  3  through  5  are  replaced 
by  the  following  two  steps: 


Step  3.5  -  The  Poisson  A  Ml1  PF  is  generated  from  the  AMU 


raraneter  in  lino  1,  col  mn  4  and  stored  in  lino  9. 

Step  4.5  -  The  AMI’  PF  (lino  9 )  and  the  AMP  PF  (line  P) 
are  conhined  as  mutually  exclusive  events  to  pet  a  resultant 
TOTMl’  Pr.  The  probabilities  of  all  possible  combinations 
are  added  together  to  get  the  probability  of  a  certain 
number  of  demands  on  the  particular  part  for  any  reason 
(such  as  missing  subassemblies  or  the  part  being  broken 
itself).  These  combinational  probabilities  are  calculated 
as  follows: 

P(0  demands)  *  P(0  AMU)  x  P(0  A’.’P) 

P(1  demand)  -  P  ( 1  AMU)  x  P(0  AWP)  +  P(0  AMU)  x  P(1  Al.’P) 


Using  the  sane  data  and  parts  hierarchy  as  in  the 
APPPvOX  example,  the  calculations  were  repeated  using  the 
EXACT  method.  Once  again,  the  example  starts  with  type  I 

parts  A,  B  ,  and  C.  For  part  A  there  are  no  subassemblies  so 

all  the  subassembly  BO  CDFs  are  as  follows: 

Step  1 

x  0  12  3  4 

P(('B0^x)  1.000  1  .000  1.000  1.000  1.000  .  .  . 

and  multiplying  by  columns  yields: 

x  0  12  3  4 

P(#AWP^x)  1.000  1.000  1.000  1.000  1.000  .  .  . 


Step  2,  converting  the  AVP  CDFs  to  a  PF: 


X 

0 

1 

2 

3 

4 

P(#AWP-x) 

1 .000 

o  .non 

0.000 

0.000 

0.000 

Step  3.5, 

generating,  a 

Poisson 

AMU  PF 

with  an 

AMU  = 

X 

0 

1 

2 

3 

4 

P  (  #  d  etna  nd  s 

«x)  .3679 

.  3679 

.1839 

.061  3 

.0153 

Step  4.5, 

combining  the 

AMU  PF 

and  the 

AWP  PF: 

P(TOTMU-O)  -  P(AUP  =  O)xP(AMU*0  )*=  1.0  x  .  3679  -  .  3679  . 
P(TOTMU-l)  «  P(AUP=1  )xP(AMU  =  0)  +  P(AWP  =  O)xP(AMU=*l  ) 

-  (0  x  .3679)  +  (1.0  x  .3679) 

-  .3679 


x  0 

1 

2 

3 

4 

P  ( t  d  etna  nd  s  =  x  )  .3679 

.  3679 

.1839 

.061  3 

.0153 

Step  6,  converting  this 

PF  to 

a  CDF  : 

x  0 

1 

*» 

3 

4 

P  (  #  d  etna  nd  s^x  )  .  3679 

7  3  5  8 

.9197 

.9810 

.9963  . 

Step  7,  adjusting  this 

CDF  for 

a  stock 

level 

of  0  and 

of  1  gives: 

x  0 

1 

2 

3 

4 

Pf/'HO^x)  .3679 

73  58 

.9197 

.9810 

.9963  . 

Step  9  ,  converting  it  to  a  R  O  P  F : 


X 

0 

1 

2 

3 

P(r«RO-X) 

.  3679 

.  3679 

.1839 

.061  3 

Step  9 

Variance 

*  1.0 

Mean  or 

F.  BO  »  1. 

0 

Variance 

to  mean 

ratio 

=  1.0 

Step  10,  the  RO  CDF  In  step  7  is  transferred 
H  for  use  in  its  calculations. 

For  part  R  the  step  results  are  as  follows: 


Step  1 

X 

0 

1 

0 

4m 

3 

P(  ^  AVJ  P-^x  ) 

1.000 

1.000 

1 . 000 

1 .000 

Step  2 

X 

0 

1 

2 

3 

p(#av;p«x) 

1 .000 

0.000 

0.000 

0.000 

Step  3.5,  generating 

a  Poisson 

AMU  PF 

for  AM 

X 

0 

1 

2 

3 

P ( # d  enand s * x ) 

.1353 

.2707 

.  2707 

.1805 

Step  4.5 

x 

0 

1 

2 

3 

P(#demands*x) 

.1353 

.  2707 

.  2707 

.1805 

4 

.0153  . 


to  parent 


4 

1.000  .  . 


4 

0.000  .  . 

-  2: 

4 

.0902  .  . 


4 

.0902  .  . 


46 


part 


&&&&&&&& 


Stop  6 


x  0 

1 

2 

3 

4  t 

P ( f d ena nd s^x  )  .1353 

.4060 

.  6767 

.  8  57  1 

.9474  .  .  . 

Step  7,  adjusting  for 

a  stock 

level  of 

1  and 

a  OPP.A  of  2: 

x  0 

1 

2 

3 

4 

P("R0^x)  .4060 

.  8  57  1 

.9834 

.9989 

.9999  .  .  . 

Step  8,  converting  to 

a  PF  : 

x  0 

1 

2 

3 

4 

P(*R0»x)  .4060 

.451  1 

.1263 

.0155 

.0011  .  .  . 

Step  9  Variance  *  .  5442 

Mean  or  ERO  =  .7546 

Variance  to  mean  ratio  *  .7211 

Step  10,  transfer  the  P  0  CDF  in  step  7  to  part  11. 

For  part  C 

Step  1 

x  0  12  3  4 

P(l»R0^x)  1.000  1.000  1.000  1.000  1.000  .  .  . 

Step  2 

x  0  12  3  4 

P()»B0-x)  1.000  0.000  0.000  0.000  0.000  .  .  . 


Step  3.5,  renoratinc  a  Poisson  A ".l'  PF  with  A‘T  =  1: 


X 

0 

1 

3 

4 

P ( # d en andsax) 

.  3679 

.  367  9 

.18  39 

.0613 

.0153  .  .  . 

Step  4.5,  TOTMU  PF : 

X 

0 

1 

2 

3 

4 

P( t A ena nd  s*x  ) 

.3679 

.  3679 

.1839 

.  061  3 

.0153  .  .  . 

Step  6,  TOTMU 

CDF  : 

X 

0 

1 

2 

3 

4 

P ( # denands^x ) 

.  3679 

.  7  3  5  8 

.9197 

.9810 

.9963  .  .  . 

Step  7,  adjusting  the 

CDF  for 

stock 

level  of 

O  and  a  0PHA  of 

2: 

X 

0 

1 

2 

3 

4 

P(^BO^x) 

.  3679 

.9197 

.9963 

.9999 

1.000  .  .  . 

Step  8,  BO  PF  : 

X 

0 

1 

2 

3 

4 

P(#BO-x) 

.  3679 

.518 

.0766 

.0036 

.0001  .  .  . 

Step  9  Variance  =  .3790 

Mean  or  EBO  =  .7162 

Variance  to  mean  ratio  *  .5292 

Step  10,  transfer  the  BO  CDF  from  step  7  to  part  H. 

Nov  H,  the  type  II  part,  is  considered. 


8 


Stop  1 


X 

0 

1 

2 

3 

4 

F(/fA  RO^x) 

.  3670 

.  73  5? 

.91^7 

.9810 

.  9963 

p  ( *  p  p,  o  <:  x  ) 

.4060 

.  P  5  7  1 

.0834 

.9989 

.9999 

P(#C  RO^x) 

.  367  9 

.9197 

.9963 

.9999 

1  .000 

P(//)l  AWP^x) 

.0549 

.  5800 

.9012 

.9799 

.  9963 

Step  2,  AW P  PF 

: 

X 

0 

1 

2 

3 

4 

P(/<AWP-x) 

.0549 

.  5  25  1 

.3212 

.0787 

.0164 

Step  3.5,  generating 

a  Poisson 

AMU  PF 

with  AMU 

=  1: 

X 

0 

1 

2 

3 

4 

P(  //de mands“x) 

.3679 

.  3679 

.1839 

.0613 

.0153 

Step  4.5,  TOTMU  PF  : 

X 

0 

l 

2 

3 

4 

P  (  #  d  enia  nd  s“x  ) 

.0202 

.2134 

.3214 

.  2471 

.127  1 

Step  6,  TOTMU 

CDF  : 

X 

0 

1 

2 

3 

4 

P(^demands^x) 

.0202 

.  2  336 

.  5  550 

.8020 

.  9292 

Step  7,  adjusted  CDF 

for  stock 

level 

of  2  and 

a  OPHA 

X 

0 

1 

-> 

3 

4 

P(0F O^x) 

.  5550 

.  802  1 

.  9292 

.  9787 

.  994  5 

Step  8,  BO  PF  : 


x  0  12  3  4 

P(/'BO-=x)  .5550  .2471  .1271  .0496  .0157  .  .  . 

Step  9  Variance  ■  1.059 

Mean  or  FRO  ■  .7422 
.Variance  to  mean  ratio  =  1.427 

Step  10,  transfer  to  parent  of  II,  if  there  was  one.  There 
is  no  parent  so  the  desired  result  is  the  ERO  of  part  H. 

In  this  simple  parts  hierarchy,  both  the  EB0  and  the 
variance  of  the  part  H  BO  PF  were  increased  by  using  the 
Poisson  approximation.  The  EP.O  was  increased  16%  (from 
.7422  to  .8612)  and  the  variance  was  increased  36.6%  (from 
1.059  to  1.4466). 

Sample  computer  output  for  the  above  example  is 
included  in  appendix  D  for  both  programs. 

Verification  o f  Computer  Models 

To  insure  APTROX  and  EXACT  function  correctly  the  per¬ 
formances  of  each  individual  subprogram  and  the  total  pro¬ 
gram  were  evaluated. 

The  R'ACT  function  was  checked  to  insure  it  produced  the 
proper  factorials  over  the  entire  ranee  of  positive 


integers.  The  results  of  the  C E N P F  subroutine  were  checked 
against  tabulated  Poisson  values  for  several  MUs.  The 
conversion  subroutines  PF2CDF  and  CDF2PF,  the  adjustment 
subroutine  SHIFTL,  and  the  combining  subroutine  COMPYN  were 
also  checked  for  proper  operation. 

The  overall  operation  of  the  programs  was  checked  by 
comparing  the  computer  output  with  hand  calculated  data  for 
the  simple  sample  problem  in  this  chapter.  The  programs 
were  also  run  for  the  examples  given  in  various  AFLC/XRS 
reports  (Ref  9,15).  All  results  agreed  to  at  least  four 
decimal  places. 


Ill  Findings 


The  analysis  was  conducted  in  two  parts.  The  first 
part  was  to  determine  what  actually  causes  the  errors  in  the 
Poisson  approximation  and  to  verify  the  conclusions  drawn  by 
constructing  parts  hierarchies  and  data  sets  which  produced 
no  error.  The  second  part  of  the  study  examined  the  sensi¬ 
tivity  of  the  errors  to  different  paramenters. 

Causes  of  Error 

To  determine  the  exact  source  of  the  approximation 
error  many  different  data  sets  were  run  through  both  the 
APPROX  and  EXACT  computer  programs.  Most  cases  exhibited 
some  error  but  in  some  cases  no  error  was  observed.  An  in 
depth  analysis  resulted  in  the  following  conclusions. 

The  Poisson  approximation  initially  preserves  the 
correct  mean  but  can  distort  the  variance  of  the  TOTMU  PF . 
The  distorted  variance  can  subsequently  induce  errors  for 
the  following  two  reasons: 

1.  The  approximation  may  assume  a  different  variance  for 
the  TOTMU  PF  than  it  actually  has,  and  when  shifted  for 
stock  level  and  QPHA,  it  distorts  the  true  mean  (this 
error  does  not  occur  when  the  part  has  a  stock  level  of 
zero  and  a  QPHA  of  one  because  no  shifting  of  the  TOTMU 
PF  takes  place).  This  error  is  referred  to  as  type  1 
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error  in  this  report. 


2  .  Also  as  a  result  of  the  Poisson  approximation,  the  P. 0 
PFs  may  have  more  or  less  variance  than  they  should 
(whether  shifted  of  not).  When  the  BO  CDFs  for  SAs  are 
transferred  to  their  HA  and  multiplied  column  bv  column 
to  get  an  AMP  CDF  for  the  HA,  the  mean  is  again  dis¬ 
torted  hy  the  incorrect  variance  even  if  the  means  or 
E  ROs  were  correct  to  start  with.  (This  error  does  not 
occur  when  there  is  only  one  SA  per  HA,  because  no  mul¬ 
tiplying  of  BO  CDFs  by  columns  tabes  place.  In  this 
case,  the  HA  AWP  CDF  is  simply  the  SA  RO  CDF.)  This 
error  is  referred  to  as  type  2  error  is  this  report. 

It  should  also  be  pointed  out  that  no  error  is  ever 
induced  for  tie  parts  at  the  very  tips  of  the  hierarchial 
tree.  Because  all  parts  are  assumed  to  have  a  Poisson  dis¬ 
tributed  AMU  and  the  tip  parts  have  no  SAs,  both  APPROX  and 
EXACT  carry  the  same  distributions  up  the  tree  to  the  first 
II  As.  It  is  then  at  this  level  that  the  error  is  first 
introduced  . 

To  show  that  the  Poisson  approximation  preserves  the 
correct  mean  but  can  distort  the  variance,  consider  the  fol¬ 
lowing  two  PFs  as  an  example. 


AN  AMU  PF 


x  0  12 

P ( X*x )  .  A  .3  .1 

for  this  PF  mean  =  .5  and  variance  *.45 

AN  AW P  PF 

x  0  12  3 

P ( X=x )  .1  .3  .4  .2 

for  this  PF  mean  *1.7  and  variance  *.81 

Using  the  APPROX  method  of  combining  PFs  the  means  are 
added  to  get  a  mean  of  2.2  from  which  the  following  Poisson 
PF  is  generated. 

RESULTANT  APPROX  TOTMU  PF 
x  0  1  2  3  4  5 

P(X-x)  .11  .24  .27  .20  .11  .05 

for  this  PF  mean  *2.2  and  variance  *2.2 


54 


^sV 


..S'-, 


L'sin"  r  h c  I' V A C T  npthn^ 


P(  »0  ) 

m 

.6 

X 

.  1 

m 

.06 

P(X-1  ) 

• 

.  3 

X 

.  1 

+ 

.6  x  .3 

m 

.  2  1 

P ( X-2  ) 

m 

.  l 

X 

.  1 

+ 

.  3  x  .3 

+ 

.6  x  .4  -  .34 

P  (  X  -  3  ) 

m 

.  l 

X 

.3 

+ 

.3  x  .4 

+ 

.6  x  .2  -  .27 

P  (  X  -  4  ) 

m 

.  l 

X 

.4 

+ 

.3  x  .2 

m 

.10 

P  (  X-  5  ) 

m 

.  l 

X 

.2 

m 

.02 

which  yields  the  following  PF : 

RESULTANT  EXACT  TOTMU  PF 
x  0  1  2  3  A  5 

P(X»x)  .06  .21  .34  .27  .10  .02 


for  this  PF  mean  -2.2  and  variance  -1.26 


This  example  illustrates  that  the  approximation 
the  same  mean  but  distorts  the  variance. 

Working  a  similar  example  in  general  terns  using 
following  two  PFs : 


PF  1 

x  0  12 

P(X-x)  ARC 


mean 


B  +  2C 


P(X*x)  a  h  c 

mean  ■  h  +  2c 

Using  the  APPROX  method  of  combining  PF s  the  means  are 
added  to  get  a  mean  of  B  +  b  +  20  +  2c. 

Using  the  EXACT  method 
P(X-O)  -  Aa 
P(X-l)  «  Ab  +  Ba 
P(X~2)  -  Ac  +  Bb  +  Ca 

P(X-3)  -  Be  +  Cb 
P(X«4)  -  Cc 

from  which  the  following  expression  for  the  mean  is  obtained 
using  expected  value: 

Mean  -  Ab  +  Ba  +  2Ac  +  2Bb  +  2Ca  +  3Bc  +  3Ch  +  4Cc 

Since 

A  +  B  +  C  »  1  and  a  +  b  +  c  *  1 

it  follows  that 

A  *  1  -  B  -  C  and  a  *  1  -  b  -  c 

Substituting  these  values  into  the  above  equation  for  the 


mean  vields 


The  expressions  for  the  ATP P. n v  anf*  nxACT  neans  are  the 
sane,  d emons t r a t i n g  that  t^e  mean  is  initially  preserved 
with  the  approximation. 

Type  2.  Error.  To  illustrate  the  type  1  error  caused  by 
the  shifting  of  a  TOTMU  PF  for  stock  level  and  OPHA,  con¬ 
sider  the  demand  distributions  given  below  for  a  certain 
part.  Roth  PFs  have  the  same  mean,  but  the  first  PF  has  a 
larger  variance,  such  as  could  have  been  induced  by  the 
approximation  . 

APPROX  TOTMU  PF 
x  0  12  3  4 

P(X-x)  .2  .2  .3  .2  .1 

mean  *  1.8  variance  ■  1.56 

EXACT  TOTMU  PF 
x  0  12 

P(X-x)  .1  .3  .35 

mean  *  1.8  variance  ■  1.06 

If  the  distributions  are  adjusted  for  a  stock  level  of  one. 


3  4 

.2  .05 


APPP.OX  BO  PF 


x  0  12  3  4 

P ( X  =  x )  .4  .3  .2  .1  0 

with  mean  *  1  and  variance  =  1 

EXACT  BO  PF 

x  0  12  3  4 

P(X-x)  .4  .35  .2  .05  0 

with  mean  *  .9  and  variance  *  .79 

In  this  case  shifting  the  distribution  for  a  stock  level  of 
one  results  in  the  approximation  giving  a  high  mean  and 
preserving  the  inflated  variance. 

If  the  original  distributions  are  adjusted  for  a  QPHA 
of  two,  the  BO  PFs  are: 

APPROX  BO  PF 

x  0  12  3  4 

P(X-x)  .2  .5  .3  0  0 

with  mean  *  1.1  and  variance  ■  .49 


EXACT  B  0  PF 


x  0  12  3  4 

P(X-x)  .1  .65  .25  0  0 

with  mean  *  1.15  and  variance  *  .33 

In  this  case,  shifting  the  distribution  for  a  OPHA  of  tw 
results  in  the  approximation  giving  a  low  mean  and  preserv 
ing  the  inflated  variance. 

If  the  original  distributions  are  adjusted  for  a  OPH 
of  two  and  a  stock  level  of  one,  the  BO  PFs  are: 

APPROX  BO  PF 

x  0  12  3  4 

P(X-x)  .4  .5  .1  0  0 

with  mean  *  .7  and  variance  ■  .41 

EXACT  BO  PF 

x  0  12  3  4 

P(X-x)  .4  .55  .05  0  0 

with  mean  »  .65  and  variance  ■  .33 


In  this  final  case,  shifting  the  distribution  results  in  the 


approxination  giving  a  high  mean  and  preserving  the  inflated 
variance . 


This  example  shows  that  the  distorted  variance  caused 
hy  the  Poisson  approximation  induces  an  error  in  the  mean, 
and  that  this  error  can  be  understated  as  well  as  over¬ 
stated. 


According  to  the  explanation  of  type  1  error,  any 
single-level  parts  hierarchy  having  an  HA  with  a  QPHA  of  1 
and  a  stock  level  of  zero  should  have  no  error  induced, 
regardless  of  the  AMU  of  the  HA  or  the  AML’s,  stock  levels, 
or  OPIIAs  of  the  subassemblies. 


This  hypothesis  was  tested  with  several  different  data 
sets,  and  there  was  never  any  error  observed.  However,  when 
any  stock  level  was  included  for  the  TIA,  or  the  OPilA  was 
changed  from  one,  errors  resulted.  Percent  error  was  used 
to  measure  both  mean  and  variance  error  and  was  calculated 
uslne  the  following  formula: 


APPROX  results  -  EXACT  results 


Percent  error  * 


EXACT  resul ts 

The  data  sets  used  and  the  results  are  included 


in  appendix 


E  . 


Type  Error.  To  illustrate  the  type  2  error  caused  by 


transferring 


PO  CDFs  for  f. As  to  their  respective  HA,  con- 


siller  the  demand  distributions  riven  below  for  a  certain 
oart.  Roth  PFs  have  the  same  no  an,  but  the  first  bas  a 
larger  variance,  such  as  could  have  been  induced  bv  the 
a  nnr  ox ira  a  t i o  n  . 


APPROX  BO  PFs  for  2  S As 


0  12  3 


P  (  X  **  x  )  .2  .4  .2  .2 


with  mean  ■  1.4,  and  variance  *  1.04 


EXACT  BO  PFs  for  2  SAs 


P ( X*x )  .1  .5  .3  .1 


with  mean  =  1.4,  and  variance  m  .64 


Combining  the  two  BO  PFs  in  each  case  by 

1.  converting  the  PFs  to  CDFs 

2.  squaring  the  cumulative  probabilities  in  each  column 

3.  converting  the  resultant  CDF  back  to  a  PF 
yields  the  following  AWP  PFs  for  the  HA: 


APPROX  AWP  PF  for  HA 


P(X-x)  .04  .32  .28  .36 


x>.v  v.v.vy.v..‘.v v .v . 


with  mean 


1.96,  and  variance 


84 


EXACT  AWP  PF  for  HA 

x  0  12  3 

P(X=x)  .01  .35  .45  .19 

with  mean  =  1.82,  and  variance  =  .55 

This  example  shows  that  the  inflated  variance  initially 

caused  by  the  approximation  can  distort  the  mean  in  later 
calcul ations . 

To  test  the  type  2  error  hypothesis,  multiple  indenture 
level  hierarchies  with  one  SA  per  HA  on  all  but  the  lowest 
level  were  run.  Again  a  QPHA  of  one  and  a  stock  level  of 

zero  for  all  parts  were  used  to  prevent  any  type  1  error. 

As  expected,  no  error  was  recorded  until  more  than  one  SA 

was  associated  with  an  HA,  above  the  lowest  level.  The 
results  of  this  test  are  included  in  appendix  F. 

It  can  therefore  be  concluded  that  whenever  an  HA  has 

nonzero  stock  level  or  a  QPHA  not  equal  to  one,  or  the 

hierarchy  has  more  than  one  SA  per  HA  above  the  lowest  level 
in  each  branch,  error  may  be  induced. 


62 


Sensitivity  Analysis 


Sensitivity  analysis  was  performed  to  study  the  effect 
of  the  following  parameters  on  the  Poisson  approximation 
error: 

-  depth  of  indent  tire  of  the  parts  hierarchy 

-  SA  to  HA  ratio  for  higher  levels 

-  stock  level  of  each  part 

-  QPIIA  of  each  nart 

To  check  the  effect  of  increasing  depth  of  indenture, 
four  different  parts  hierarchies  were  constructed,  ranging 
from  a  " s  i  ng le- 1  eve  1 ”  model  to  a  "four-level"  model.  To 
insure  that  the  change  in  the  error  was  due  only  to  increas¬ 
ing  depth  of  indenture,  the  same  basic  part  data  was  used  in 
each  additional  level.  The  hierarchies,  part  data  sets,  and 
results  are  included  in  appendix  C.  A  summary  of  the 
results  is  shown  below. 

ERROR  OF  THE  MEAN 
THE  POISSON  APPROXIMATION 

Depthof  Setl  Set  2 

Indenture  7  error  %  error 

single-level  0  0 

two-level  1.4  -7.6 

three-level  3.3  -15.2 

four-level  5.0  -22.3 

The  results  suggest  that  the  magnitude  of  the  error  in  the 
mean  for  the  Poisson  a ppr ox ima t i o n  increases  with  increasing 
depth  of  indenture.  This  is  expected  since  deeper  indenture 


increases  the  nnnher  of  opportunities  for  type  1  and  tvpe  2 
error  to  occur. 

To  evaluate  the  effect  of  SA  to  HA  ratio,  four  dif¬ 
ferent  parts  hierarchies  were  constructed.  These  hierar¬ 
chies  all  had  two  levels  of  indenture  and  had  constant  SA  to 
HA  ratios  ranging  from  two  to  five  (for  example,  for  a  con¬ 
stant  SA  to  HA  ratio  of  three,  each  HA  in  the  hierarchy 
would  have  three  SAs).  These  parts  hierarchies  were  all  run 
for  four  different  data  sets.  The  hierarchies,  part  data 
sets,  and  results  are  included  in  appendix  C.  A  summary  of 
the  results  is  shown  below. 


ERROR  OF  THE  MEAN  US  I  MO 

Tin:  poissou  approximation 


SA  to  HA 
Ratio 

Set  1 
%  error 

Set  2 
%  error 

Set  3 
terror 

Set  4 
terror 

2 

1.5 

-7.6 

2.8 

-6.7 

3 

3.2 

-6.9 

6.2 

-8.7 

4 

5.0 

-5.6 

8.6 

-10.0 

3 

6.8 

-3.P 

10.5 

-12.0 

The 

results  suggest 

that  the 

error  in 

the 

mean 

for  the  Po  is- 

son 

approxima  tion 

can  move  in 

either  direction 

with  increas- 

i  ng 

SA  to  HA  ratio 

• 

This  is 

probably 

due 

to  the  fact  that 

a  s 

the  SA  to  HA 

ratio  increases  type 

2 

error 

tends  to  don- 

inatc .  Depending 

on 

whether 

tvpe  2  error 

i  s 

positively  or 

negatively  inclined,  the  error  moves  in  the  plus  or  minus 
direction. 


To  study  the  effect  of  different  combinations  of  stock 


Level 


QPHA 

2 

3 

5 . 5 

5.0 

A. 3 

4.2 

1  .  6 

1  .  3 

.  6 

- .  2"  -  ^ 

.  1 

The  lack  of  monotonicity  in  the  above  results  makes  it 
difficult  to  arrive  at  specific  conclusions.  However,  there 
are  general  trends  present  in  all  three  cases  which  may  shed 
some  light  on  the  problem.  In  all  cases  there  seems  to  be: 

a.  A  general  decrease  in  positive  error  (not  magni¬ 
tude  of  error)  with  increasing  stock  level 

b.  A  general  increase  in  positive  error  (not  magni¬ 
tude  of  error)  with  increasing  QPHA  (in  the  low 
stock  level/high  QPHA  sectors  this  does  not  hold). 

c.  A  general  cancelling  of  errors  between  the  two 
factors  as  evidenced  by  the  negative  errors  in  the 
high  stock  and  low  QPHA  sectors,  and  by  the 
minimum  errors  occurring  roughly  along  the  indi¬ 
cated  diagonals. 

d.  The  errors  of  hierarchy  two  are  generally  less 
than  the  errors  of  the  other  two  hierarchies. 


.v.v.v.v. v  ■.■v/.-v  .  ■. 


IV  Conclusions  and  Recommendations 


The  WARS  function  of  building  optimal  WRM  kits  is 
extremely  important.  Stocking  insufficient  quantities  could 
severely  degrade  the  ability  to  successfully  fight  a  war, 
while  stocking  excess  quantities  or  unnecessary  parts  is 
very  expensive  in  holding  and  opportunity  costs.  This  is 
graphically  portrayed  by  the  fact  that  a  WRSK  for  a  single 
24  aircraft  A-10  squadron  costs  in  excess  of  $8  million  (Ref 
17).  The  ability  of  WARS  to  build  optimal  kits  is  dependent 
on  its  ASSESS  subroutine  being  able  to  accurately  evaluate 
the  performance  of  different  parts  packages  in  the  marginal 
analysis  process.  Since  ASSESS  uses  the  Poisson  approxima¬ 
tion,  understanding  and  evaluating  its  effect  is  of  utmost 
importance . 

Evaluating  the  effect  of  this  approximation  is 
extremely  difficult  because  of  the  infinite  number  of  possi¬ 
ble  hierarchy,  MU,  stock  level,  and  QPHA  combinations,  and 
the  interaction  of  the  above  parameters.  Several  specific 
hierarchies  and  data  sets  can  easily  be  evaluated,  but  using 
induction  to  extrapolate  this  limited  set  of  results  to  the 
general  case  is  dangerous.  Another  uncertainty  is  the  real¬ 
ism  of  the  samples  evaluated.  All  parts  hierarchies  used 
were  fabricated  since  the  part  indenture  relationships 
required  are  not  currently  available  for  operational  air- 


craft.  Determining  the  effect  of  an  approximation  without 
actual  data  is  a  formidable  problem,  nevertheless,  from  this 
study  the  following  conclusions  were  drawn: 


Conclusions 

(Note:  Mean  error  means  the  error  of  the  mean  not  average 
error . ) 


1.  Mean  error  usually  is  induced  when  the  Poisson 
approximation  is  used. 


2.  The  approximation  initially  preserves  the  correct 
mean  but  can  distort  the  variance  of  the  part's 
demand  probability  function. 


3.  The  distorted  variance  of  a  part's  demand  proba¬ 
bility  function  can  induce  errors  in  the  mean  in 
subsequent  computations. 


4.  The  errors  in  the  mean  can  be  caused  by  two  dis¬ 
tinct  processes: 


Type  1  -  distortion  of  the  mean  when  a  PF  is 
shifted  for  stock  level  or  QPHA  because  the 
approximation  may  assume  a  different  variance 
for  the  TOTMU  PF  than  actual,  (this  error 
does  not  occur  when  the  part  has  a  stock 
level  of  zero  and  a  QPHA  of  one  because  no 
shifting  of  the  TOTMU  PF  takes  place). 

Type  2  -  distortion  of  the  mean  when  the  BO 
CDFs  for  SAs  are  transferred  to  their  HA  and 
multiplied  column  by  column  to  get  an  AWP  CDF 
for  the  HA,  because  the  Poisson  approximation 
distorts  the  variance  jf  the  BO  PFs  (This 
error  does  not  occur  when  there  is  only  one 


SA  per  HA,  because  no  multiplying  of  BO  CDFs 
by  columns  takes  place.  In  this  case,  the  HA 
AWP  CDF  is  simply  the  SA  BO  CDF.) 


5.  As  a  result  of  type  1  error,  an  error  may  be 
Induced  whenever  stock  level  is  more  than  zero  or 
QPHA  is  greater  than  one  for  any  parent  part. 

6.  As  a  result  of  type  2  error,  an  error  may  also  be 
induced  whenever  the  parts  hierarchy  is  such  that 
an  HA  has  more  than  one  SA  (excluding  the  bottom 
level)  . 

7.  If  neither  of  the  above  conditions  in  conclusion  5 
or  conclusion  6  are  present,  no  error  will  result. 

8.  Both  mean  and  variance  error  vary  in  sign  and  mag¬ 
nitude  from  case  to  case. 

9.  Mean  error  seems  to  increase  in  magnitude  with 
increasing  depth  of  indenture. 

10.  Mean  error  seems  to  increase  or  decrease  monotoni- 
cally  with  increasing  SA  to  HA  ratio,  depending  on 
the  sign  of  the  original  type  2  error. 

11.  The  interaction  of  stock  level  and  QPHA  makes  it 
difficult  to  arrive  at  specific  conclusions.  How¬ 
ever,  there  are  general  trends  present  in  all 
three  cases  which  may  shed  some  light  on  the  prob¬ 
lem.  In  all  cases  there  eeems  to  be: 


70 


a.  A  general  decrease  in  positive  error  (not 

magnitude  of  error)  with  increasing  stock 
level 

b.  A  general  increase  in  positive  error  (not 

magnitude  of  error)  with  increasing  QPHA  (in 
the  low  stock  level/high  QPHA  sectors  this 
does  not  hold ) . 

c.  A  general  cancelling  of  errors  between  the 
two  factors  as  evidenced  by  the  negative 
errors  in  the  high  stock  and  low  QPHA  sec¬ 
tors,  and  by  the  minimum  errors  occurring 
roughly  along  the  indicated  diagonals. 

d.  The  errors  of  hierarchy  two  are  generally 
less  than  the  errors  of  the  other  two  hierar¬ 
chies.  This  suggests  that  minimally  inden¬ 
tured  hierarchies  induce  less  error  than 
heavily  indentured  hierarchies. 

The  error  created  may  not  be  critical  to 

WARS/ASSESS  because: 

a.  There  may  be  some  cancelling  of  errors,  both 
between  type  1  and  type  2  error  as  well  as 
between  the  errors  caused  by  the  other 
approximations  used  in  ASSESS. 


b.  Most  actual  aircraft  parts  hierarchies  have  a 
small  percentage  of  their  parts  indentured  to 
higher  assemblies  (Ref  8,18).  Because  of 
this  characteristic,  the  approximation  errors 
are  diluted  somewhat  and  tend  to  be  smaller. 

c.  ASSESS  is  only  used  in  the  marginal  analysis 
process  to  evaluate  the  Increase  in  perfor¬ 
mance  due  to  adding  one  part  to  a  kit  rela¬ 
tive  to  another  different  part.  Since  both 
kits  evaluated  would  be  the  same  kit  except 
for  a  few  different  parts,  they  would  most 
likely  suffer  from  similar  errors  and  the 
methodology  would  probably  still  be  able  to 
pick  the  part,  which  if  added,  would  increase 
performance  most  per  dollar. 

Recommendations 

Although  this  study  produced  only  a  limited  number  of 
specific  conclusions,  it  is  hoped  that  the  data  computed  and 
the  general  conclusions  presented  aid  in  the  further  under¬ 
standing  of  this  Important  approximation.  To  build  upon 
these  results,  the  following  areas  for  further  study  are 
suggested : 

1.  In  the  future,  when  indentured  information  is 
available  for  operational  aircraft,  it  would  be 


valuable  to  run  actual  hierarchies  and  parts  data 
through  the  APPROX  and  EXACT  programs  and  analyze 
the  results. 

Similar  studies  should  be  undertaken  analyzing  the 
effects  of  the  other  approximations  used  in 
WARS/ ASSESS . 

The  combined  effects  of  the  three  approximations 
should  be  examined  in  detail. 

Since  WARS  will  have  the  capability  to  use  nega¬ 
tive  binomial  distributions  for  each  part's  AMU 
distribution  instead  of  Poisson,  the  study  should 
be  repeated  for  this  distribution.  The  Poisson 
distribution  is  a  special  case  of  the  negative 
binomial  with  the  variance  to  mean  ratio  equal  to 
one,  so  sensitivity  to  the  variance  to  mean  ratio 
should  also  be  evaluated. 
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Appendix  A 


Glossary  (Ref  10) 

AMU  -  A  variable  denoting  the  basic  mean  pipeline  quantity 
for  an  item,  excluding  AWP.  Also  called  "actual  mu." 

Assessment  -  In  this  report,  assessment  is  the  calculation 
of  the  expected  NMCS,  total  expected  backorders,  and 
probability  distribution  of  NMCS  based  on  a  specific 
stock  level. 

Availabibli ty  -  The  probability  that  a  randomly  selected 
item  is  not  backordered. 

AWP  -  A  variable  denoting  the  expected  number  of  items 
"awaiting  parts.”  This  quantity  is  computed  by  the  WARS 
assessment  routine. 

Backorder(  BO)  -  An  unfilled  demand  or  shortage  at  Mae 
level . 

BLSS  kit  -  Stock  of  “Base  Level  Self-sufficiency  Spares" 
designed  to  meet  the  additional  requirements  for  spare 
parts,  above  normal  peacetime  requirements,  when  a  war 
starts. 

Canmble  -  Interchangeable  with  canniballzable . 

Ca nnibal izable  item  -  An  item  which  can  be  borrowed  from  an 
aircraft  that  already  has  a  backorder  to  prevent  a 
backorder  on  another  aircraft.  The  extra  maintenance 
time  required  to  cannibalize  a  part  is  not  considered, 
so  an  item  which  takes  an  excessive  amount  of  time  to 
cannibalize  should  be  considered  non-cannibal izable  . 
Canniballzable  is  synonymous  with  cannable. 

CDF  -  Cumulative  Distribution  Function 

Combined  probabilities  -  This  refers  to  calculating  the  pro¬ 
bability  distribution  of  the  sum  of  probability  distri¬ 
butions  for  two  mutually  exclusive  random  variables. 
For  example,  denoting  two  random  variables  by  A  and  B, 
and  the  probability  of  x  demands  for  A  by  p(A  *  x),  the 
probability  distribution  of  the  sum  of  demands  for  A 
and  B  is  calculated  as  follows: 
p(A+B-0)  -  p(A-0)p(B-0) 
p(A+B-l)  -  p(A-0)p(B-l)  +  p(A-l)p(B-0) 


Cumulative  probability  -  The  probability  of  0  through  x 
demands,  since  only  discrete,  non-  negative  distribu¬ 
tions  are  considered. 

DFHP  -  Daily  Flying  Hour  Program.  The  planned  number  of 
hours  a  unit  will  fly  each  day  of  a  contingency. 

Demand  -  Demand  is  interchangeable  with  "pipeline  quantity" 
and  failure. 

Depth  of  Indenture  -  refers  to  the  number  of  levels  in  a 
parts  hierarchy.  The  more  levels  in  a  parts  hierarchy, 
the  greater  the  depth  of  indenture. 

EBO  -  Expected  number  of  backorders.  If  the  item  under  con¬ 
sideration  is  the  aircraft  itself,  EBO  is  equivalent  to 
E [ NMCS ] . 

E(NMCS]  -  Expected  number  of  NMCS  aircraft.  Also  called 
E[NMC]  in  this  report. 

Failure  -  Failure  is  interchangeable  with  demand  and  "pipe¬ 
line  quantity." 

HA  -  Higher  Assembly. 

Hierarchy  diagram  -  A  diagram  showing  the  indenture  rela¬ 
tionships  for  all  positions  on  the  aircraft.  For  each 
item,  this  Indenture  relationship  lists  all  higher 
assemblies  that  must  be  taken  apart  to  reach  this  item, 
and  all  lower  indenture  items,  or  parts  that  can  only 
be  reached  by  taking  this  item  apart. 

Hierarchy  position  -  The  location  of  a  part  on  the  hierarchy 
diagram.  A  hierarchy  position  describes  the  location 
of  only  one  stock  number,  but  there  may  be  more  than 
one  unit  of  the  stock  number  in  the  hierarchy  poslton 
if  its  quantity  per  aircraft  is  larger  than  one.  A 
stock  number  will  have  more  than  one  hierarchy  position 
(and  will  be  called  a  non-unique  item)  if  it  is 
atttached  to  more  than  one  parent  item  on  the  aircraft. 

Higher  assembly  (HA)  -  An  item  that  contains  parts  that  can 
only  be  reached  by  taking  apart  this  item.  Sometimes 
called  a  "parent"  with  respect  to  immediately  lower 
parts  . 

Imperfectly  cannibalizable  part  (MCP)  -  An  item  that  is  con¬ 
sidered  cannbal i zable  but  does  not  have  the  same  quan¬ 
tity  per  higher  assembly  on  all  higher  assemblies. 
Also  referred  to  as  "myopically  cannibalizable”. 


Indenture  -  The  concept  of  determining  which  items  must  be 
taken  apart  to  reach  a  given  part.  An  item  is  top 
indenture  if  it  can  be  removed  directly  from  the  air¬ 
craft,  regardless  of  where  it  is  attached  to  the  air¬ 
craft.  A  lower  Indenture  part  can  be  reached  only  by 
taking  apart  the  item  to  which  it  is  attached. 

Individual  probability  -  the  probability  of  exactly  x 
backorders  (or  other  demands). 

Kit  -  A  stock  of  spare  parts  built  for  use  at  a  single  squa¬ 
dron,  regardless  of  whether  the  stock  consists  of  a 
WRSK,  a  BLSS,  or  normal  peacetime  supplies. 

Lower  indentured  item  -  A  part  whcih  can  be  reached  only  by 
taking  apart  an  item  to  which  this  part  is  attached. 

LRU  -  Line  Replaceable  Unit.  An  item  which  can  be  removed 
and  replaced  on  the  flightline.  It  contains  one  or 
more  lower  Indenture  items. 

MA  -  Marginal  Analysis. 

MCP  -  Imperfectly  cannibalizable  ,  myopically  cannibal izable  , 
or  Imperfectly  cannable  part. 

MU  -  pipeline  quantity. 

Myopic  cannibalization  -  Same  as  Imperfect  cannibalization. 

NCP  -  Non-Cannibal i zable  Part. 

NMC  -  Not  Mission  Capable,  used  interchangeably  with  NMCS. 

NMCS  -  Not  Mission  Capable  due  to  Supply.  -  An  aircraft 
missing  parts  is  called  NMCS,  but  a  spare  part  can  only 
be  AWP  if  it  has  lower  items  backordered.  Previously 
called  NORS  (Not  Operationally  Ready  due  to  Supply). 

Non-cannibal izable  part  (NCP)  -  An  item  that  cannot  be 
removed  in  a  serviceable  condition  from  one  aircraft  to 
prevent  a  backorder  on  another  aircraft.  While  it  may 
be  physically  possible  to  remove  such  an  item,  the  item 
may  be  designated  NC  because  the  item  would  be  damaged 
when  it  is  removed  or  time  to  remove  it  is  considered 
excessive. 

Non-unique  -  A  stock  number  that  is  used  in  more  than  one 
hierarchy  position  on  an  aircraft. 

NSN  -  National  Stock  Number. 
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Parent  item  -  The  next  higher  assembly  directly  above  an 
item  in  an  indenture  relationship. 

PCP  -  Perfectly  Cannibal izable  Part 

Perfectly  c annibal i za bl e  part  (PCP)  -  An  item  which  is  con¬ 
sidered  cannibal i zab le ,  and  which  has  the  same  quantity 
per  higher  assembly  on  all  higher  assemblies  containing 
the  item. 

PF  -  Probability  function. 

Pipeline  quantity  (MU)  -  The  number  of  units  of  an  item  that 
have  been  removed  from  aircraft  and  are  not  yet  in  a 
serviceable  condition  at  the  base  ready  to  be  replaced 
on  aircraft.  This  includes  parts  that  are  awaiting 
base  maintenance,  in  base  maintenance,  awaiting  parts 
at  the  base,  awaiting  or  in  transportation  to  the 
repair  depot,  awaiting  maintenance  at  the  depot,  in 
maintenance  at  the  depot,  and  awaiting  or  in  transpor¬ 
tation  to  the  base.  Demand,  failure,  or  backorder  may 
be  used  as  a  substitute  for  expected  pipeline  quantity. 

Probability  function  (PF)  -  The  set  of  the  individual  proba¬ 
bility  of  each  possible  outcome  for  a  discrete  random 
variable . 

QPA  -  Quantity  Per  Aircraft.  This  is  the  number  of  each 
item  on  one  aircraft. 

QPHA  -  Quantity  Per  Higher  Assembly.  This  is  the  number  of 
each  item  on  its  respective  higher  assembly. 

SA  -  Subassembly. 

SDO  -  Stock  Due  Out.  The  total  number  of  parts  on 

backorder . 

SRU  -  Shop  Replaceable  Unit.  An  item  which  can  be  reached 
only  by  removing  the  LRU  containing  the  item  from  the 
aircraft  and  disassembling  the  LRU. 

Stock  Level  -  the  initial  number  of  each  part  available  to 

replace  failed  parts.  For  WRSK  and  BLSS  it  is  the 

quantity  included  in  the  kit.  For  peacetime  operation 

it  is  the  quantity  in  base  supply. 

Stock  number  -  Interchangeable  with  NSN. 

Subassembly  (SA)  -  A  part  indentured  to  a  higher  assembly. 

TOTMU  -  A  variable  denoting  the  total  expected  pipeline 
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quantity  for  an  item.  It  Is  equal  to  AMU  +  AWP. 

Type  (of  an  Item)  -  The  type  of  an  Item  refers  to  its  inden¬ 
ture  relationship  as  follows: 

4  ■  An  LRU  or  engine  attached  directly  to  the  aircraft. 

3  ■  An  LRU  or  module  attached  to  a  type  4  item. 

2  ■  An  LRU  attached  to  a  module. 

1  *  An  SRU  attached  to  any  LRU. 

UE  -  Unit  equipment. 

Unit  equipment  (UE)  -  The  number  of  aircraft  assigned  to  a 
unit. 

Unique  Item  -  A  stock  number  used  in  only  one  hierarchy 
position. 

WARS  -  Wartime  Assessment  and  Requirements  System. 

WRSK  -  War  Readiness  Spares  Kit.  A  kit  of  spare  parts  to  be 
taken  with  a  unit  in  wartime,  designed  to  support  that 
unit  for  a  specified  period  without  outside  maintenance 
support . 


Ap  pe  nd  lx  B 


The  Assessment  Rout ine  (Ref  £) 

This  section  discusses  the  assessment  routine  as  it  is 
used  in  WARS  ASSESS.  First,  some  terms  need  to  be  defined: 

Item  types: 

Type  I  -  an  SRU  assigned  to  an  LRU. 

Type  II  -  an  LRU  assigned  to  a  module. 

Type  III  -  an  LRU  or  module  assigned  to  an 

engine. 

Type  IV  -  an  LRU  or  engine  assigned  to  an 
aircraf t . 

Stages  one  through  four  of  the  computation  cover  Type 
I  through  Type  IV  items  respectively. 

i  ■  present  item  considered,  in  squadron  S  on  day  t. 

It  may  be  n  stages  denoted: 

i  “  item  which  will  not  be  cannibalized, 
n 

i  *  item  which  is  perfectly  cannibalized, 

pc 

i  ■  item  which  is  imperfectly  cannibalized, 
m  c 

H(i)  »  higher  assembly  for  item  i. 

J  ■  quantity  of  i  in  an  H(i). 

J  “  maximum  quantity  of  i  in  any  H(i)  ■ 

J(H(i) ,S,t J . 

K  ■  quantity  of  H(i)  considered  at  a  certain  step 


in  Che  process 


MU[(i,S,t)]  *  M  U (  i )  *  expected  number  of  assets  of 
item  i  accounted  for  at  the  end  of  day  t  in  squadron  S 
(basically  equals  the  number  of  demands  that  must  be  filled 
by  base  stock  to  prevent  backorders). 

p  (x , MU [ i , S  ,  t  ]  )  *  p(x,MU[i])  *  probability  of  x  assets 
accounted  for,  given  an  expected  number  of  MU  [  i , S , t  ]  . 

A[i,S,t]  *  A[i]  *  number  of  assets  of  item  i  con¬ 
sidered  at  the  end  of  day  t  in  squadron  S  (basically  equals 
the  number  of  serviceable  items  in  base  stock). 

N[  j,H(i),S,t]  -  N[j,H(i)]  -  number  of  H(i)  that  con¬ 
tain  j  of  item  i  at  the  end  of  day  t  in  squadron  S. 

N  l  ~>A  » H(  * )  t  S  ,  t  ]  ■  N[^j,H(i))  ■  number  of  H(i)  that  con¬ 
tain  between  j  and  J  of  item  i  at  the  end  of  day  t  in  squa¬ 
dron  S  . 


N[H(i) ,S,t]  *  N[H(i)]  *  total  number  of  higher  assem¬ 
blies  considered  in  this  run. 

This  may  include  some  higher  assemblies  without  item 
i,  or  N(0,H(i)].  Normally  the,  number  of  higher  assemblies 
considered  will  be  the  number  installed  on  all  aircraft  in 
the  squadron. 

EBO(i,S,t]  ■  E  BO [  i  ]  ■  expected  number  of  backorders  of 
item  i  in  squadron  S  at  the  end  of  day  t. 


Qfll(l)  ,  S  ,  t  ]  -  Q  [  H  ( i )  ]  -  QPH  A  of  item  i  on  H(i). 

F[K,i,S,t]  ■  F [ K , i ]  “  probability  that  K  units  or  less 
of  H(i)  are  missing  item  i  at  the  end  of  day  t  in  squadron 
S,  assuming  i  is  perfectly  c annib al i za b le  . 

P(K,H(i) ,S,t]  -  P [ K , H ( i ) ]  -  probability  that  K  units 

of  H(i)  are  missing  one  or  more  i  items  at  the  end  of  day 

pc 

t  in  squadron  S. 

q [  i  ,  S  ,  t ]  ■  q[i]  ■  expected  proportion  of  units  of  H(i) 
that  are  not  missng  any  of  item  i  at  the  end  of  day  t  in 
squadron  S,  assuming  i  is  non-cannibalizable. 

C [ K  ,  i , S , t ]  -  C[K,i]  -  probability  that  K  units  of  H(i) 
are  missing  item  i  at  the  end  of  day  t  in  squadron  S,  assum¬ 
ing  that  i  is  Imperfectly  cannibalizable ,  and  where  canni¬ 
balization  consolidates  the  failures  onto  the  smallest 
number  of  H( i )  . 

M[H(i),S,t]  ■  M[H(i)]  *  expected  number  of  H(i)  miss¬ 
ing  one  or  more  of  item  i  at  the  end  of  day  t  in  squadron  S, 
assuming  i  is  imperfectly  cannibalizable. 

g[l,S,t]  ■  gli]  *  expected  proportion  of  units  of  11  (i) 
that  are  not  missing  any  of  item  i  at  the  end  of  day  t  in 
squadron  S,  assuming  1  is  Imperfectly  cannibalizable. 

b[H(l),S,t]  »  b I H ( ± )  )  *  the  probability  that  a  random 
unit  of  H(l)  is  missing  any  i  or  i  items  at  the  end  of 


day  t  in  squadron  S 


B l K , H ( i) , S , t ]  *  B[K,H(i)]  ■  an  approximate  (binomial) 

probability  that  K  units  of  II ( i )  are  missing  one  or  more  i 

or  i  at  the  end  of  day  t  in  squadron  S. 
me 

4>  [  K  ,H(  i  )  ,  S  ,  t  ]  ■  $[K,H(i)]  ■  An  approximate  probability 
that  K  units  of  H(i)  are  down  (or  NMCS)  at  the  end  of  day  t 
in  squadron  S. 

E[H(i),S,t]  *  E [ H ( i ) ]  ■  the  expected  number  of  units 
of  H ( i )  down  (or  NMCS)  at  the  end  of  day  t  in  squadron  S. 

The  procedure  to  be  used  in  the  example  below  is  as 
follows:  Steps  a  through  d  are  performed  once,  at  the 

beginning  of  the  run,  for  all  items.  Each  stage  of  the  com¬ 
putation  starts  with  step  e. 

Step  a.  Partition  all  items  into  the  categories  i  , 

pc 

i  ,  and  i 
n  me 

Step  b.  For  each  item,  input  A[i]. 

Step  c.  For  each  item,  determine  N[j,H(i))  for  all 
appropriate  J. 

Step  d.  For  each  item,  compute  N[H(i)]  ■ ^  N [ j  ,  H ( 1 ) ] . 

j 

It  is  best  to  include  j  ■  0  for  those  H(i)  that  could 
contain  i  but  do  not.  Then,  N[H(i))  will  be  the  sane  for 


all  lower  Indentured  Items  i  with  in  the  same  H(i). 


Step  e.  For  each  item  in  this  stage,  compute  MU [ i ] . 
Actually,  an  "initial  MU"  will  be  calculated  for  each  item 
at  the  beginning  of  the  run.  Ths  "initial  MU"  is  the 
expected  pipeline  quantity  for  the  item,  assuming  no  shor¬ 
tages  of  lower  indentured  items  and  no  stock  of  this  item. 

In  stages  2,3,  and  4,  the  "final  MU",  or  MU[i]  as  used 
here,  will  equal  the  "initial  MU“  plus  E[H(i)J,  which  is 
calculated  for  item  i  in  step  p  of  the  previous  stage.  This 
means  that  MU[i]  is  equal  to  the  expected  pipeline  quantity 
including  shortages  of  lower  indentured  items.  For  example, 
supppose  i  in  stage  2  and  H(i)  in  stage  1  refer  to  item  LI, 
and  the  “initial  MU"  for  LI  is  0.5.  This  means  that  an 
average  of  0.5  units  of  LI  have  been  removed  from  the  air¬ 
craft  and  have  not  yet  been  returned  to  a  serviceable  condi¬ 
tion  at  the  base,  excluding  those  units  of  LI  that  are 
repaired  but  missing  lower  indentured  parts,  because  of 
shortages.  Assume  that  these  lower  indentured  parts  are 
SRUs  SI  and  S2  with  MU [ S 1 ]  -  0.5  and  MU(S2]  -  1.0.  Also 
assume  that  in  stage  1,  E[Ll]  is  found  to  be  1.3,  meaning 
that  the  failures  of  an  average  of  0.5  +  1.0  •  1.5  units  of 
SI  and  S2  can  be  consolidated  onto  1.3  units  of  LI.  Assum¬ 
ing  that  the  failures  of  LI,  excluding  those  awaiting  lower 
indentured  parts,  are  not  consolidated  with  the  units  of  LI 
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that  are  awaiting  lower  indentured  parts,  M U [ L 1 ]  ■  0.5  +  1.3 

-  1.8. 


Step  f.  For  each  item  in  this  stage,  compute  expected 
backorders,  or 

CO 

EBOfi]  -  ^  <x  “  A[i])  P<*,MU[i])  <1> 

x-A( i] 

Actually,  the  expected  backorders  do  not  need  to  be 
computed  for  c a nnibal i zab 1 e  items  unless  that  data  is  needed 
for  the  objective  function  or  for  other  evaluation  purposes. 
A  more  efficient  formula  may  be  used  for  computer  implemen¬ 
tation: 

CO 

EBO [ i ]  -  ^  (1  -  p( x , MU [ i ]  ) 

x-A{i] 

In  this  formula,  the  upper  limit  may  be  replaced  by 
A[i]  plus  the  number  of  parts  installed  on  all  the  aircraft 
in  the  squadron. 


Step  g.  For  each  i  item,  compute: 

pc 


A  [  I  ]  +  K  •  Q  [  H  (  i  )  ] 

FlK,il  “  2  p(  x  ,  MU  {  i  ]  ) 

x  *  0 

for  K  -  0,1,2, ...  ,(N[H(i)  -  N[0,H(i)]) 


<  2  > 


Step  h.  For  each  i  item,  compute: 


q{i]  -  >  n  1 J , n ( i )  1  1 1  -  eboqj _ \  j  < 3 > 

^  j^N  [  H  (  i )  ]  ^  2  j  •  N[j,H(i))  J 

Step  i.  For  each  i  item  compute  C ( K ( i ) ]  for  all  K  1 

me 

0,1,2,...,  (N  [H(  i)  ]  -  N[0,ll(i)J): 


A  [  i  ]  +  J  •  K 


for  K<  N [ J , H (  i  ) ]  ,  C { K 


p ( x  ,  MU [ i ]  ) 


for  N [ J , H ( i ) ]  <  K  4  N[>J-1,  H(i>], 

A [ i ]  +  J  •  N { J , H (  i  )  ]  +  (J~l )(K-N[J , H ( i ) ]  ) 


C  [  K  ,  i  ]  - 


p( x , MU { i ]  ) 


for  N[^J-l,H(i)]  <  K  N<  N[>J-2,H(i)  ]  , 

A[i]+J«  N[J,H<i) ]+<J-l)N{J-l,H(i))  +  (J~2)(K-N[J,H(i)]) 


C[K,i] 


■2 


p( x , MU [ i ] ) 


for  N[>2,H(i))  <  K  4  -Nf>l,H(i)J, 

A[  i]+J*  N[  J,H(i)  ]+.  .  .+2  N[2  ,H(i)  ]  +  l  ( K-N  [  >,2  , H (  i  )  ]  ) 


-  •  2 


p ( x , MU [ i]  ) 


In  other  words,  this  means  to  calculate  the  cumulative 
probability  of  less  than  or  equal  to  K  units  of  H(i)  down 


for  all  relevant  K,  assuming  that  cannibalization  consoli¬ 
dates  the  failures  onto  the  smallest  possible  number  of 
H(i).  This  happens  by  cannibalizing  in  order  from  the  H(i) 
containing  the  most  i  to  the  least  i. 

Then  , 


M[H(i)  ] 


N  [  H  (  i  )  -  N  [  0  ,  H  (  i  )  ] 


2 


K(C[K,i] 


C [ K- 1  ,  i  ]  ) 


Compute  : 

g [ i 1  "  1  -  M[H(i)  /  N [ H (  i  )  ]  <4> 

Step  j.  For  all  i  items  in  each  H(i),  compute: 

n 

n  '  i  <5> 

n 

i 

n 

The  q  values  are  from  formula  <3>. 

Step  k.  For  all  i  items  in  each  H(i)  compute: 

m  c 

n  st1  i  <6> 

me 

i 

me 

The  g  values  are  from  formula  <4>. 


Step  1.  For  each  H(i)  compute: 


b i h ( i ) ]  « 1  -  n  <ii 1 1  •  n  *! 1  i 


n  '  me 

i  i 

n  me 

Step  m.  For  each  H(i),  compute  binomial  probabilities 
for  all  K  -  0,1,2, ...,NfH(i)l: 

K  x  N  [  H ( i ) ] -x 


N[H(i)] ! 


(N[H(i)]-x) !x! 


b[H(i>]  (l-b|H(i)))  <7  > 


x-0 

Step  n.  For  all  K  *  0 , 1 , 2 , . . . , N [ H ( i ) )  compute  for 

each  H( i) : 

P[K,H(i)]  -  J[  F[K,i  ]  <8> 

pc 

i 

pc 

The  F  values  are  from  formula  <2>. 

Step  o.  For  all  K  ■  0 , 1 , 2 , . . . N ( H ( i ) ]  compute  for  each 

H(i): 


0[K,H(i)]  -  (B{K,H(i))(P[K,H(i)]) 


<9  > 


Step  p.  For  each  H(i)  compute: 

N(H(i)] 

E [ H ( i )  ]  -  ^  K($(K,H(i)]  -  *[K-l,H(i)]) 

K-l 


<10> 


Formula  <10>  is  the  number  of  H(i)  expected  to  be  down 


due  to  lack  of  lower  indentured  items.  Actually,  in  an 
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operational  implementation,  a  more  efficient  formula  would 
be  used: 


E  J  H( i)  ] 


NfH(i) ]-l 

N[H(i)]  -  ^  0[K,H(i)] 

K»0 


At  the  end  of  stage  4,  the  value  of  E[H(1>]  is  the 
expected  number  of  NMCS  aircraft,  since  H(i)  is  now  the  air¬ 
craft.  If  desired,  the  total  number  of  expected  backorders 
i  s : 


EBO[i].  All  results  apply  to  squadron  S  at  time  t 


Appendix  C 


Program  Listings 


Note:  The  programs  are  written  in  FORTRAN  77  and  can  be  run  for  any 
data  set  by  simply  changing  the  DATA  statement  and  the  NPARTS  parameter 
throughout  the  programs.  A  list  of  computer  variables  used  follows  the 
two  program  listings. 


PROGRAM  APPROX 


******************************************************************* 

*  THIS  PROGRAM  CALCULATES  THE  APPROXIMATE  E[NMC1  AND  INTERMEDIATE  * 

*  DEMAND  AND  BACKORDER  DISTRIBUTIONS  FOR  A  SPECIFIED  PARTS  * 

*  HIERARCHY  AND  STOCK  LEVEL  * 

******************************************************************* 

I NTEGER  P , L , C , TYPE , BLOCK , NPARTS , DIS  TL 
CHARACTER*40  MSG(1:20) 

PARAMETER  (NPARTS-3 .DISTL-2 7, BLOCK-6 *(NPARTS+1 )*(DISTL+1 ) ) 
REAL  A(0:NPARTS,1:14,0:DISTL) 

COMMON  A 

DATA  ( ( A(P , 1 , C) ,C-1 , 7 ) ,P-0 .NPARTS ) / 

+  5, 1,1,0. 0,0,0, 

+  4, 3, 1,1, 2,0,1, 

+  3,0,1, 1, 0,1,1, 

+  3, 0,2, 2, 1,1, 2, 

+  3,0, 2, 1,0, 1,3/ 

+  ( ((A(P,L,C),P-0, NPARTS), L-2, 7), C-O.DISTL)/ 

+  BL0CK*1 .0/ 

+  MSC/'0/TYPE/#SA/QPHA/MU/SL/HA#/PUP' ,'SA  1  BO  CDF', 

+  'SA  2  BO  CDF ','SA  3  BO  CDF’.'SA  4  BO  CDF’.'SA  5  BO  CDF’, 

+  ' AWP  CDF'.'AWP  PF' ,'VAR/AWP/VAR  TO  MEAN  RATIO/TOTMU' , 

+  'TOTMU  PF', 'TOTMU  CDF', 

+  'BO  CDF', 'BO  PF' ,' VAR/MEAN/VAR  TO  MEAN  FOR  BO, TOTMU, &  AWP', 
+  ******  ERROR  ****** , ' IMPROPER  PART  TYPE  IN  DATA’, 

+  ’NEGATIVE  VALUE  IN  DATA',’>  5  SA  PER  PARENT  PART', 

+  ’HA*  >  NPARTS', 'PUP  >  PARENT ”S  #SA'/ 


*  ECHO  PRINT  DATA  AND  ERROR  CHECK 
PRINT*, 'INPUT  DATA' 

PRINT* 

PRINT*, 'PART  TYPE  #SA  QPHA  MU  SL  HA#  PUP' 
DO  10  P-0, NPARTS 

WRITE(*,5)P,(A(P, 1 ,C),C-1,7) 

5  FORMAT(I5,7F6.1) 

IF  (A(P,1,1).GT.5.0  .OR.  A( P, l , 1) .LT. 1 .0)  THEN 
PRINT*,MSG(I5) 


PRINT*, MSG(16) 

END  IF 

IF  (A(P,1,2).GT.5.0)  THEN 
PRINT*, MSG( 15) 

PRINT*, MSG(18) 

END  IF 

IF  (A(P, 1 ,6) .GT.NPARTS)  THEN 
PRINT*, MSG(15) 

PRINT*, MSG(19) 

END  IF 

IF  (A(P, 1 ,7) .GT. A( A(P, 1 ,6) ,1,2))  THEN 
PRINT*, MSG( 15) 

PRINT*, MSG( 20) 

END  IF 
DO  20  C-1,7 

IF  (A(P,1,C).LT.0.0)  THEN 
PRINT*, MSG(15) 

PRINT*,MSG(17) 

END  IF 

20  CONTINUE 

10  CONTINUE 
PRINT* 

STEP  THROUGH  INDENTURE  LEVELS  FROM  LOWEST  TO  HIGHEST  LEVEL 
DO  40  TYPE-1,5 

FOR  EACH  PART,  IF  IT  IS  THE  TYPE  BEING  CONSIDERED,  PERFORM  CALCULATIONS 
DO  50  P-0 , NPARTS 
IF  (A(P, 1 , 1 ) .EQ.TYPE)  THEN 

COMPUTES  AWP  CDF  FOR  PART 

DO  60  C-O.DISTL 
DO  70  L-2,6 

A(P,7,C)-A(P,7,C)*A(P,L,C) 

70  CONTINUE 

60  CONTINUE 

CONVERT  AWP  CDF  TO  AWP  PF 

CALL  CDF2PF(P, 7 ,8) 

COMPUTE  AWP  AND  VARIANCE;  ADD  AWP  TO  AMU 
IF  PART  IS  THE  AIRCRAFT  THEN  CALCULATE  AWP  AND 
AWP  VARIANCE  AND  GO  TO  THE  END  OF  THE  LOOP 
IF  (TYPE.EQ.5)  THEN 
CALL  CALEB0(P,8, 14, 1 ) 

CALL  CALVAR(P,8, 14,0) 

GO  TO  40 
ELSE 

CALL  CALEBO(P,8,9,l) 

A(P,14,15)-A(P,9,1) 

IF  (A(P,9,1) .GT.O.O)  THEN 
CALL  C ALVAR( P , 8 , 9 , 0 ) 


*************************************************** A*** ******** ******** 


FUNCTION  FACT(Y) 

*  THIS  FUNCTION  CALCULATES  THE  FACTORIAL  OF  A  GIVEN  ARCUMENT 
INTEGER  Y,K 
REAL  FACT 
IF  (Y.LT,0)  THEN 

PRINT*,' ERROR-TRIED  TO  TAKE  FACTORIAL  OF  A  NEGITIVE  NUMBER* 

END  IF 
FACT-1 

DO  100  K-l  ,Y 
FACT«FACT*K 
100  CONTINUE 
END 

*********************************************************************** 


SUBROUTINE  PF2CDF(P,LF,LT) 

*  THIS  SUBROUTINE  CONVERTS  A  PF  TO  A  CDF 
INTEGER  P,LF,LT,K,DISTL,NPARTS 
PARAMETER  (NPARTS-3 , DISTL-27 ) 

REAL  A(0: NPARTS, 1: 14,0:DISTL) 

COMMON  A 

A(P,LT,0)-A(PtLF ,0) 

DO  110  K-l.DISTL 

A(P,LT,K)-A(P,LF,K)+A(P,LT,K-l) 

110  CONTINUE 
END 

*********************************************************************** 


SUBROUTINE  CDF2PF(P,LF,LT) 

THIS  SUBROUTINE  CONVERTS  A  CDF  TO  A  PF 


s 

INTEGER  P,LF,LT,K,DISTL, NPARTS 

p 

PARAMETER  (NPARTS-3, DISTL-27) 

REAL  A( 0: NPARTS, 1: 14, 0:DISTL) 

COMMON  A 

m 

A(P,LT,0)-A(P,LF,0) 

i 

DO  120  K-l.DISTL 

A(P, LT,K)-A(P,LF,K)-A(P#LF,K-1 ) 
120  CONTINUE 
END 


********************************************************************** 


SUBROUTINE  GENPF(P,MU,LT) 

*  THIS  SUBROUTINE  GENERATES  A  POISSON  PROBABILITY  FUNCTION 
INTEGER  P,LT,Y,DISTL, NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  MU, FACT 

REAL  A(0: NPARTS, I : 14,0:DISTL) 


COMMON  A 
DO  130  Y-0, DISTL 

A(P,LT, Y )«MU*  *Y*EXP (-MU ) / F ACT ( Y ) 

130  CONTINUE 
END 

*********************************** A*************** ******* ************ 


SUBROUT INE  SHIFTL(PF ,LF , PT , LT , SL, QP ) 

*  THIS  SUBROUTINE  ADJUSTS  A  DEMAND  CDF  FOR  STOCK  LEVEL  AND 

*  QPHA  TO  GET  A  BACKORDER  CDF 

INTEGER  PF , LF , PT , LT , SL, QP , K , DISTL , NPARTS 
PARAMETER  (NPARTS-3, DISTL-27) 

REAL  A(0: NPARTS, 1 : 14,0: DISTL) 

COMMON  A 

DO  140  K-0,INT((DISTL-SL)/QP) 

A( PT , LT, K)-A( PF , LF , SL+K*QP ) 

140  CONTINUE 

DO  145  K«INT((DISTL-SL)/QP)+1, DISTL 
A(PT,LT,K)»1.0 
145  CONTINUE 
END 

********************************************************************** 


SUBROUTINE  CALEBO( P , LF , LT , CT ) 

*  THIS  SUBROUTINE  CALCULATES  EXPECTED  BACKORDER  FROM  A  PF 
INTEGER  P,LF,LT,CT,K, DISTL, NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  EBO 

REAL  A(0:NPARTS,1:14,0:DISTL) 

COMMON  A 
EBO-O.O 

DO  150  K-0, DISTL 

EBO-EBO+K*A(P,LF,n 
150  CONTINUE 
A(P,LT,CT)-EB0 
END 

********************************************************************** 


SUBROUTINE  COMBYN(P , LF 1 , LF2 , LT ) 

*  THIS  SUBROUTINE  COMBINES  TWO  MUTUALLY  EXCLUSIVE  PFS  TO  GET 

*  A  RESULTANT  PF 

INTEGER  N,C,P,LF1,LF2,LT,N, DISTL, NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  A(0: NPARTS, 1:14,0: DISTL) 

COMMON  A 
DO  160  C-0, DISTL 
A(P,LT,C)»0 .0 
DO  160  N-0,C 

A(P,LT,C)-A(P,LT,C)+A(P,LF1,N)*A(P,LF2,C-N) 


160 


CONTINUE 


END 

********************************************************************* 


SUBROUTINE  CALVAR(P,LF,LT,CT) 

*  THIS  SUBROUTINE  CALCULATES  THE  VARIANCE  AND 

*  VARIANCE  TO  MEAN  RATIO  OF  A  BACKORDER  PF 

INTEGER  P,LF,LT,CT,K,DISTL,NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  SUM 

REAL  A(0:NPARTS,1:14,0:DISTL) 

COMMON  A 
SUM-0.0 

DO  170  K-0,DISTL 

SUM«SUM+( (K-A(P,LT, CT+1 ))**2)*A(P,LF ,K) 

170  CONTINUE 
A(P,LT,CT)-SUM 

A(P,LT,CT+2)-A(P,LT,CT)/A(P,LT,CT+l) 

END 

********************************************************************* 


PROGRAM  EXACT 


******************************************************************* 

*  THIS  PROGRAM  CALCULATES  THE  EXACT  E[NMC]  AND  INTERMEDIATE  * 

*  DEMAND  AND  BACKORDER  DISTRIBUTIONS  FOR  A  SPECIFIED  PARTS  * 

*  HIERARCHY  AND  STOCK  LEVEL  * 

******************************************************************* 


I NTEGE R  P , L , C , TYPE , BLOCK , NPARTS , DI STL 
CHARACTER*40  MSG(1:20) 

PARAMETER  (NPARTS-3 .DISTL-27 , BLOCK-6* (NPARTS+1 )*(DISTL+1 ) ) 
REAL  A(0:NPARTS,1:14,0:DISTL) 

COMMON  A 

DATA  ((A(P,1,C),C-1, 7), P-0, NPARTS)/ 

+  5, 1,1, 0,0,0, 0, 

+  4, 3, 1,1, 2, 0,1, 

+  3, 0,1,1, 0,1,1, 

+  3, 0,2, 2, 1,1, 2, 

+  3, 0,2, 1,0, 1,3/ 

+  (((A(P,L,C), P-0, NPARTS), L-2, 7), C-O.DISTL)/ 

+  BL0CK*1 .0/ 

+  MSG/’0/TYPE/#SA/QPHA/MU/SL/HA#/PUP* ,'SA  1  BO  CDF’, 

+  'SA  2  BO  CDF’.'SA  3  BO  CDF ’,'SA  4  BO  CDF','SA  5  BO  CDF’, 

+  'AWP  CDF', 'AWP  PF'.'AMU  PF’.'TOTMU  PF’.’TOTMU  CDF', 

+  *  BO  CDF’,  'BO  PFV  VAR/ME  AN  /VAR  TO  MEAN  FOR  B0,T0TMU,6i  AWP', 

+  ******  ERROR  *****',' IMPROPER  PART  TYPE  IN  DATA’, 

+  ’NEGATIVE  VALUE  IN  DATA’,’>  5  SA  PER  PARENT  PART’, 

+  ’HA#  >  NPARTS’, 'PUP  >  PARENT ’ ’ S  #SA’/ 


*  ECHO  PRINT  DATA  AND  ERROR  CHECK 
PRINT*, ’INPUT  DATA’ 

PRINT* 

PRINT*, ’PART  TYPE  #SA  QPHA  MU  SL  HA#  PUP’ 
DO  10  P-0, NPARTS 

WRITE(*,5)P,(A(P,1 ,C) ,C-1 ,7) 

5  FORMAT(I5,7F6.1) 

IF  (A(P,1,1).GT.5.0  .OR.  A(P, 1 , l ) .LT. 1 .0)  THEN 
PRINT*, MSG( 15) 

PRINT* ,MSG(16) 

END  IF 

IF  (A(P,1,2) .GT.5.0)  THEN 
PRINT*, MSG( 15) 

PRINT*, MSG( 18) 

END  IF 

IF  (A(P,1, 6). GT. NPARTS)  THEN 
PRINT*,MSG(15) 

PRINT*, MSG(19) 

END  IF 


IF  (A(P,1,7).GT.A(A(P,1,6),1,2))  THEN 
PRINT*,MSC(15) 

PRINT*, MSG( 20) 

END  IF 
DO  20  C-1,7 

IF  <A(P,1,C).LT.0.0)  THEN 
PRINT*, MSG( 15) 

PRINT*,MSG(17) 

END  IF 

20  CONTINUE 

10  CONTINUE 
PRINT* 

*  GENERATE  AMU  DISTRIBUTION  FOR  EACH  PART  FROM  PIPELINE  QUANTITY 

DO  30  P-l.NPARTS 

CALL  GENPF(P,A(P, 1 ,4) ,9) 

30  CONTINUE 

*  STEP  THROUGH  INDENTURE  LEVELS  FROM  LOWEST  TO  HIGHEST  LEVEL 

DO  40  TYPE-1,5 

*  FOR  EACH  PART,  IF  IT  IS  THE  TYPE  BEING  CONSIDERED,  PERFORM  CALCULATIONS 

DO  50  P-O.NPARTS 
IF  ( A( P , 1 , 1 ) . EQ . TYPE )  THEN 

*  COMPUTES  AWP  CDF  FOR  PART 

DO  60  C-O.DISTL 
DO  70  L-2,6 

A(P,7,C)-A(P,7,C)*A(P,L,C) 

70  CONTINUE 

60  CONTINUE 

*  CONVERT  AWP  CDF  TO  AWP  PF 

CALL  CDF2PF(P,7,8) 

*  CALCULATE  MEAN  AWP  AND  AWP  VARIANCE 

*  IF  PART  IS  THE  AIRCRAFT  THEN  GO  TO  THE  END  OF  THE  LOOP 

IF  (TYPE.EQ.5)  THEN 
CALL  CALEBO( P , 8 , 14 , 1 ) 

CALL  CALVAR(P,8,14,0) 

GO  TO  40 
ELSE 

CALL  CALEBO (P, 8, 14, 15) 

IF  (A(P,14,15).GT.0.0)  THEN 
CALL  CALVAR(P,8, 14,14) 

END  IF 
END  IF 

*  COMBINE  AWP  PF  AND  AMU  PF  TO  GET  TOTMU  PF 

*  COMPUTE  MEAN  AND  VARIANCE  OF  TOTMU  PF 

CALL  COMBYN(P,8,9, 10) 

CALL  CALEBO(P, 10, 14,8) 

IF  (A(P, 14,8) .GT. 0.0)  THEN 


CALL  CALVAR(P, 10, 14 , 7 ) 

END  IF 

*  CONVERT  TOTMU  PF  TO  TOTMU  CDF 

CALL  PF2CDF(P, 10, 11 ) 

*  ADJUST  TOTMU  CDF  FOR  STOCK  LEVEL  AND  QPHA  TO  GET  BACKORDER  CDF 

CALL  SHIFTL(P,11,P,12,INT(A(P,1,5)), 

+  INT(A(P,1,3))) 

*  CONVERT  BACKORDER  CDF  TO  BACKORDER  PF 

CALL  CDF2PF(P,12,13) 

*  CALCULATE  EXPECTED  VALUE  OF  BACKORDER  PF 

CALL  C ALEBO( P, 13, 14,1) 

*  CALCULATE  VARIANCE  AND  VARIANCE  TO  MEAN  RATIO  OF  BACKORDER  PF 

CALL  CALVAR( P, 13, 14,0) 

*  TRANSFER  BACKORDER  CDF  TO  PARENT  RECORD 

CALL  SHIFTL(P,12,INT(A(P,1,6)), 

+  1+INT(A(P,1,7)),0,1) 

END  IF 

50  CONTINUE 

40  CONTINUE 

*  PRINT  E[NMC]  AND  PFS 

PRINT* 

PRINT*, ’EXACT  RESULTS’ 

DO  80  P-0, NP ARTS 
PRINT* 

PRINT*, 'PART  NO.  ',P,'  **********', 

+  •  *  *************** 

DO  90  L-2,14 
PRINT* 

PRINT*, MSG(L) 

PRINT  95,(A(P,L,C) ,C-0,DISTL) 

90  CONTINUE 

95  FORMAT(7F9.5) 

80  CONTINUE 
END 

**************** A********** A* AAA**** A A ****** AAA****** A***************** 

FUNCTION  FACT(Y) 

*  THIS  FUNCTION  CALCULATES  THE  FACTORIAL  OF  A  GIVEN  ARGUMENT 

INTEGER  Y,K 
REAL  FACT 
IF  (Y.LT.O)  THEN 

PRINT*, 'ERROR-TRIED  TO  TAKE  FACTORIAL  OF  A  NEGITIVE  NUMBER' 

END  IF 


DO  100  K-1,Y 
FACT«FACT*K 
100  CONTINUE 
END 

*********************************************************************** 

SUBROUTINE  PF2CDF(P,LF,LT) 

*  THIS  SUBROUTINE  CONVERTS  A  PF  TO  A  CDF 

INTEGER  P,LF,LT,K, DISTL, NPARTS 
PARAMETER  (NPARTS«3,DISTL-27) 

REAL  A(0:NPARTS, 1 : 14,0: DISTL) 

COMMON  A 

A(P,LT,0)«A(P,LF,0) 

DO  110  K«l, DISTL 

A(P,LT,K)«A(P,LF,K)+A(P,LT,K-1) 

110  CONTINUE 
END 

*********************************************************************** 

SUBROUTINE  CDF2PF(P,LF,LT) 

*  THIS  SUBROUTINE  CONVERTS  A  CDF  TO  A  PF 

INTEGER  P ,  LF ,  LT ,  K, DISTL , NPARTS 
PARAMETER  (NPARTS-3.DISTL-27) 

REAL  A(0:NPARTS,1:14,0:DISTL) 

COMMON  A 

A(P,LT,0)-A(P,LF,0) 

DO  120  K-l, DISTL 

A(P,LT,K)-A(P,LF,K)-A(P,LF,K-1) 

120  CONTINUE 
END 

********************************************************************** 
SUBROUTINE  GENPF(P,MU,LT) 

*  THIS  SUBROUTINE  GENERATES  A  POISSON  PROBABILITY  FUNCTION 

INTEGER  P,LT,Y, DISTL, NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  MU, FACT 

REAL  A(0 : NPARTS, 1: 14,0: DISTL) 

COMMON  A 
DO  130  Y-0, DISTL 

A(P,LT,Y)-MU**Y*EXP(-MU)/FACT(Y) 

130  CONTINUE 
END 

********************************************************************** 

SUBROUT INE  SHI FTL ( PF , LF , PT , LT , SL , QP ) 

*  THIS  SUBROUTINE  ADJUSTS  A  DEMAND  CDF  FOR  STOCK  LEVEL  AND 

*  QPHA  TO  GET  A  BACKORDER  CDF 


INTEGER  PF , LF ,PT, LT, SL,QP, K, DISTL, NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  A(0: NPARTS, 1 : 14,0: DISTL) 

COMMON  A 

DO  140  K-0,INT((DISTL-SL)/QP) 

A(PT,LT,K)-A(PF,LF,SL+K*QP) 

140  CONTINUE 

DO  145  K«INT((DISTL-SL)/QP)+1, DISTL 
A(PT,LT,K)«1.0 
145  CONTINUE 
END 

********************************************************************* 

SUBROUTINE  CALEBO(P ,LF ,LT,CT) 

*  THIS  SUBROUTINE  CALCULATES  EXPECTED  BACKORDER  FROM  A  PF 

INTEGER  P,LF,LT,CT,K, DISTL, NPARTS 
PARAMETER  (NPARTS-3 .DISTL-27) 

REAL  EBO 

REAL  A(0: NPARTS, 1:14,0: DISTL) 

COMMON  A 
EBO-O.O 

DO  150  K-0, DISTL 

EBO»EBO+K*A( P , LF ,K) 

150  CONTINUE 
A(P,LT,CT)«EBO 
END 

********************************************************************** 

SUBROUTINE  C0MBYN(P,LF1 ,LF2,LT) 

*  THIS  SUBROUTINE  COMBINES  TWO  MUTUALLY  EXCLUSIVE  PFS  TO  GET 

*  A  RESULTANT  PF 

INTEGER  N,C,P,LF1,LF2,LT,N, DISTL, NPARTS 
PARAMETER  (NPARTS-3, DISTL-27) 

REAL  A(0:NPARTS,1: 14,0:DISTL) 

COMMON  A 
DO  160  C-0, DISTL 
A(P,LT,C)-0.0 
DO  160  N-0,C 

A(P,LT,C)-A(P,LT,C)+A(P,LF1,N)*A(P,LF2,C-N) 

160  CONTINUE 

END 

********************************************************************* 

SUBROUTINE  CALVAR(P,LF ,LT,CT) 

*  THIS  SUBROUTINE  CALCULATES  THE  VARIANCE  AND 

*  VARIANCE  TO  MEAN  RATIO  OF  A  BACKORDER  PF 

INTEGER  P,LF,LT,CT,K, DISTL, NPARTS 
PARAMETER  (NPARTS-3, DISTL-27) 

REAL  SUM 
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REAL  A(0:NPARTS, 1 : 14 ,0:DISTL) 

COMMON  A 
SUM-0.0 

DO  170  K-O.DISTL 

SUM-SUM+((K-A(P,LT,CT+1))**2)*A(P,LF,K) 
170  CONTINUE 
A(P,LT,CT)-SUM 

A( P , LT , CT+2 )-A( P,LT,CT)/A(P,LT, CT+1 ) 

END 


*********** A******** *********** A* *************** A* ********** A******** 


Variables  Used  in  the  Computer  Programs 


A  -  data  and  probability  array. 

APPROX  -  program  which  calculates  approximate  E[NMC). 

BLOCK  -  computed  parameter. 

BO  -  output  label  denoting  backorder. 

C  -  Indexing  variable  for  the  column  number  in  the  A  array  which  calcu¬ 
lates  the  EBO  from  a  probability  function  (PF). 

CALEBO  -  subroutine  which  calculates  the  EBO  from  a  PF. 

CALVAR  -  subroutine  which  calculates  variance  and  variance  to  mean  ratio 
of  a  backorder  PF. 

CDF2PF  -  subroutine  which  converts  a  cumulative  distribution  function 
(CDF)  to  a  PF. 

COMBYN  -  subroutine  which  combines  two  mutually  exclusive  PFs  to  get  a 
resultant  PF. 

DISTL  -  length  of  probability  distributions. 

DLPAR  -  distribution  length  paramenter. 

EBO  -  expected  backorder. 

EXACT  -  program  which  calculates  exact  E[NMC]. 

FACT  -  function  which  calculates  the  factorial  of  a  number. 

GENPF  -  subroutine  which  generates  a  Poisson  PF. 

L  -  Indexing  variable  for  the  line  number  in  the  A  array. 

LF  -  indexing  variable  indicating  the  line  number  of  array  A  from  which 
a  probability  distribution  is  taken. 

LF1  -  indexing  variable  of  array  A  indicating  the  line  from  which  the 
first  of  two  mutually  exclusive  PFs  is  taken. 

LF2  -  indexing  variable  of  array  A  indicating  the  line  from  which  the 
second  of  two  mutually  exclusive  PFs  is  taken. 

LT  -  indexing  variable  indicating  the  line  of  array  A  to  which  a  resul¬ 
tant  probability  distribution  is  assigned. 
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MSC  -  message  and  heading  array. 

MU  -  parameter  from  which  the  Poisson  distribution  is  generated  (as  an 
output  label  it  indicates  the  pipeline  quantity  for  an  item). 

NPARTS  -  the  number  of  different  parts  on  the  aircraft  hierarchy. 

P  -  indexing  variable  for  the  part  number  in  the  A  array. 

PF  -  indexing  variable  for  array  A  indicating  the  part  from  which  a  pro¬ 
bability  distribution  is  taken. 

PF2CDF  -  subroutine  which  converts  a  PF  to  a  CDF. 

PT  -  indexing  variable  for  array  A  indicating  the  part  to  which  the 
resultant  probability  distribution  is  assigned. 

PUP  -  position  of  a  subassembly  under  a  parent  item,  for  example,  if  a 
parent  item  has  two  subassemblies(SA) ,  the  first  SA  will  have  a 
PUP  of  1  and  the  second  SA  will  have  a  PUP  of  2. 

QP  -  QPHA  for  a  specific  part. 

SHIFTL  -  subroutine  which  adjusts  a  demand  CDF  for  stock  level  and  QPHA 
to  get  a  backorder  CDF. 

SL  -  stock  level  for  a  specific  part. 

SUM  -  computational  variable  used  in  calculating  variance. 

TYPE  -  part  type. 

#SA  -  output  label  indicating  the  number  of  different  subassemblies  on 
an  item. 


Appendix  D 


Sample  Results 


Note:  E[NMC]  is  the  boxed  figure  at  the  end  of  the  data  for  part  number 

0. 


Approx  Results  for  Sample  Problem 


INPUT  DATA 


PART 

TYPE 

#SA 

QPHA 

MU 

SL 

HA# 

PUP 

0 

5.0 

1.0 

1.0 

.0 

.0 

.0 

.0 

1 

4.0 

3.0 

1.0 

1.0 

2.0 

.0 

1.0 

2 

3.0 

.0 

1.0 

1.0 

.0 

1.0 

1.0 

3 

3.0 

.0 

2.0 

2.0 

1.0 

1.0 

2.0 

4 

3.0 

.0 

2.0 

1.0 

.0 

1.0 

3.0 

APPROX  RESULTS 


PART  NO. 

o  *  *  * 

*  *  *  * 

***** 

*  *  *  * 

*  *  *  *  * 

*  *  *  * 

SA  1  BO  CDF 

.54678 

.76004 

.89272 

.95875 

.98613 

.99587 

.99890 

.99973 

.99994 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  2  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

103 


/.i 


SA  5  BO  CDF 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

AWP  CDF 

.54678 

.76004 

.89272 

.95875 

.98613 

.99587 

.99890 

.99973 

.99994 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

.54678 

.21326 

.13267 

.06603 

.02739 

.00974 

.00303 

.00084 

.00021 

.00005 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/AWP/VAR  TO  MEAN  RATIO/TOTMU 


.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

BO  CDF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

BO  PF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN /VAR  TO  MEAN 
1.44662  \  .86116  jl 

FOR  BO, 
.67985 

TOTMU, & 
.00000 

AUP 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

& 


;*■> 


PART  NO. 


************************* 


SA  1  BO  CDF 

.36788  .73576 

.99999  1.00000 

1. 00000  1.00000 
1.00000  1.00000 


.91970 

1.00000 

1.00000 

1.00000 


.98101 

1.00000 

1.00000 

1.00000 


.99634 

1.00000 

1.00000 

1.00000 


.99941 

1.00000 

1.00000 

1.00000 


.99992 

1.00000 

1.00000 

1.00000 


SA  2  BO  CDF 

.40601  .85712 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 


.98344 

1.00000 

1.00000 

1.00000 


.99890 

1.00000 

1.00000 

1.00000 


.99995 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


SA  3  BO  CDF 


.36788 

1.00000 


.91970 

1.00000 


.99634 

1.00000 


.99992 

1.00000 


1.00000 

1.00000 


1.00000 

1.00000 


1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

£%- 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

V.v 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 
1.00000  1.00000 


1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


AWP  CDF 
.05495 
.99999 


.58000 

1.00000 


.90116 

1.00000 


.97985 

1.00000 


.99629 

1.00000 


.99940 

1.00000 


.99992 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

B 

AWP  PF 

.05495 

.52505 

.32116 

.07870 

.01644 

.00311 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

te: 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/AWP/VAR  TO  MEAN 

RATIO/TOTMU 

•  -  * . 

.66601 

1.48845 

.44745 

2.48845 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

*  ‘.*J 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

■ 

TOTMU  PF 
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.08304 

.20664 

.25710 

.21326 

.13267 

.06603 

.02739 

.00974 

.00303 

.00084 

.00021 

.00005 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.08304 

.28968 

.54678 

.76004 

.89272 

.95875 

.98613 

.99587 

.99890 

.99973 

.99994 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.54678 

.76004 

.89272 

.95875 

.98613 

.99587 

.99890 

.99973 

.99994 

.99999 

1.00000 

1.00000 

l.OCOOO 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.54678 

.21326 

.13267 

.06603 

.02739 

.00974 

.00303 

.00084 

.00021 

.00005 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN  FOR  BO 
1.44662  .86116  1.67985 

, TOTMU, & 
.00000 

AWP 

.00000 

.00000 

.00000 

2.48845 

2.48845 

1.00000 

.00000 

.00000 

.00000 

.00000 

.66601 

1.48845 

.44745 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO. 

2  *  *  * 

*  *  *  * 

*  *  *  *  * 

*  *  *  * 

***** 

*  *  *  * 

SA  1  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

l .00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  2  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

1.00000  1.00000 

1,00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000  1.00000  1.00000  1.00000  1.00000  1.00000  1.00000 

1.00000  1.00000  1.00000  1.00000  1.00000  1.00000  1.00000 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

l. 00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/AWP/VAR  TO  MEAN  RATIO/TOTMU 


.00000 

.00000 

.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.36788 

.73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

l. 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.36788 

.73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN  FOR  BO, TOTMU, &  AWP 


1.00000 

1.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

1.00000 

1.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO. 

3  *  *  * 

*  *  *  * 

*  *  *  *  * 

*  *  *  * 

*  *  *  *  * 

*  *  *  * 

SA  1  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  2  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  CDF 

1.00000 

l.  00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.000' I'; 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/AWP/VAR  TO  MEAN  RATIO/TOTMU 

.00000 

.00000 

.00000 

2.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 
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.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.13534 

.00344 

.27067 

.00086 

.27067 

.00019 

.18045 

.00004 

.09022 

.00001 

aBiTm  mil 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.13534 

.40601 

.67668 

.85712 

.94735 

.98344 

.99547 

.99890 

.99976 

.99995 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.40601 

.85712 

.98344 

.99890 

.99995 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.40601 

.45112 

.12631 

.01547 

.00105 

.00005 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN /VAR  TO  MEAN  FOR  BO 

,  TOTMU, & 

AWP 

.54416 

.75458 

.72114 

.00000 

.00000 

.00000 

.00000 

2.00000 

2.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO.  4  ************************* 


SA  1  BO  CDF 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  2  BO  CDF 

1.00000 

l. 00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1. 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 


1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

l. 00000 

i. 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/AWP/VAR  TO  MEAN  RATIO/TOTMU 


.00000 

.00000 

.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.36788 

.73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.36788 

.91970 

.99634 

.99992 

l. 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.36788 

.55182 

.07664 

.00358 

.00008 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

no 


.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN 

FOR  BO 

,  TOTMU , & 

AWP 

.37902 

.71617 

.52924 

.00000 

.00000 

.00000 

.00000 

1.00000 

1.00000  1 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

i 


Exact  Results  for  Sample  Problem 


INPUT  DATA 

PART  TYPE  #SA 

0  5.0  1.0 

1  4.0  3.0 

2  3.0  .0 

3  3.0  .0 

4  3.0  .0 


EXACT  RESULTS 

PART  NO.  0  *  * 

SA  1  BO  CDF 

.55499  .80204 

.99995  .99999 

1.00000  1.00000 
1.00000  1.00000 

SA  2  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  3  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  4  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  5  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 


QPHA 

MU 

SL 

1.0 

.0 

.0 

1.0 

1.0 

2.0 

1.0 

1.0 

.0 

2.0 

2.0 

1.0 

2.0 

1.0 

.0 

*  *  * 

*  *  * 

*  * 

.92915  .97872 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


HA# 

PUP 

.0 

.0 

.0 

1.0 

1.0 

1.0 

1.0 

2.0 

1.0 

3.0 

*  *  *  * 

*  *  *  * 

.99446 

.99872 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

*  *  *  * 


.99973 
1 .00000 
1.00000 
1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


.00000 

.00000 

.00000 

.00000 


AWP  CDF 
.55499 


.80204  .92915 


.97872 


.99446  .99872 


.99973 


.99995  .99999  1.00000  1.00000  1.00000  1.00000  1.00000 

1.00000  1.00000  1.00000  1.00000  1.00000  1.00000  1.00000 

1.00000  1.00000  1.00000  1.00000  1.00000  1.00000  1.00000 


AWP  PF 


.55499 

.24705 

.12711 

.04958 

.01574 

.00426 

.00101 

.00022 

.00004 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

AMU  PF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

,00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

BO  CDF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

BO  PF 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/ MEAN /VAR  TO  MEAN  FOR  BO 
1.05943  1  .74224  1 1.42734 

, TOTMU,  & 
.00000 

AWP 

. 00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

. 00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO. 

1  *  *  * 

*  *  *  * 

*  *  *  *  * 

*  *  *  * 

•k  k  *  *  * 

*  *  *  * 

SA  1  BO  CDF 

.36788  .73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000  1.00000  1.00000  1.00000  1.00000  1.00000  1.00000 


SA  2  BO  CDF 


.40601 

.85712 

.98344 

.99890 

.99995 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

.36788  .91970 

.99634 

.99992 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 . 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  CDF 

.05495 

.58000 

.90116 

.97985 

.99629 

.99940 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

.05495 

.52505 

.32116 

.07870 

.01644 

.00311 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

AMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.02021 

.21337 

.32141 

.24705 

.12711 

.04958 

.01574 

.00426 

.00101 

.00022 

.00004 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.02021 

.23358 

.55499 

.80204 

.92915 

.97872 

.99446 

.99872 

.99973 

.99995 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.55499 

.80204 

.92915 

.97872 

.99446 

.99995 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

BO  PF 

.55499 

.24705 

.12711 

.04958 

.01574 

.00022 

.00004 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN  FOR  BO 

,T0TMU,& 

AWP 

1.05943 

.74224 

1.42734 

.00000 

.00000 

1.66601 

2.48845 

.66950 

.00000 

.00000 

.66601 

1.48845 

.44745 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO. 

2  *  *  * 

*  *  ★  * 

*  *  *  *  A 

*  *  *  * 

SA  1  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

SA  2  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 


1.00000 

1.00000 

1.00000 

1.00000 

.* 

1.00000 

1.00000 

i 

" 

.99872 

.99973 

* 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

fl 

1 

.00426 

.00101 

d 

.00000 

.00000 

1 

w. 

.00000 

.00000 

.00000 

.00000 

1 

m 

4 

.00000 

.00000 

J 

.00000 

.00000 

** 

.00000 

.00000 

n 

.00000 

.00000 

fl 

i 

*  *  *  * 

*  *  *  * 

1 .00000 

1.00000 

] 

1.00000 

1.00000 

I 

1.00000 

1.00000 

1.00000 

1.00000 

4 

* 

1.00000 

1.00000 

,« 

1.00000 

1.00000 

4 

1 

1.00000 

1.00000 

1 

1.00000 

1.00000 

n 

« 

1.00000 

1.00000 

1.00000 

1.00000 

i 

1.00000 

1.00000 

■ 

1.00000 

1.00000 

* 

1 

1.00000 

1.00000 

"l 

1.00000 

1.00000 

1 

1.00000 

1.00000 

1.00000 

1.00000 

,* 

1.00000 

1.00000 

1.00000  1.00000 

1.00000  1.00000 


1 1 
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1  tl  j  1  a  1  »!■  1  fcla  1 1 


•  V-  .*•  m'- 

■  *  -  2 


1.00000  1.00000  1.00000  1.00000  1.00000  1.00000  1.00000 


AWP  CDF 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1 .00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

AMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00300 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

. 00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.36788 

.73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1. 00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.36788 

.73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN  FOR  BO 
1.00000  1.00000  1.00000 

,  TOTMU, & 
.00000 

AWP 

.00000 

.00000 

.00000 

1.00000 

1.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO.  3  ************************* 


SA  1  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  2  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  3  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  A  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 

SA  5  BO  CDF 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 


1.00000  1.00000 
1.00000  1.00000 
1.00000  1.00000 
1.00000  l. 00000 


AWP  CDF 
1.00000 
1.00000 
1.00000 
1.00000 

AWP  PF 
1.00000 
.00000 
.00000 
.00000 

AMU  PF 
.1353A 
.00344 
.00000 
.00000 

TOTMU  PF 
.13534 
.00344 
.00000 


1.00000 

1.00000 

1.00000 

1.00000 


.00000 

.00000 

.00000 

.00000 


.27067 

.00086 

.00000 

.00000 


.27067 

.00086 

.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


.00000 

.00000 

.00000 

.00000 


.27067 

.00019 

.00000 

.00000 


.27067 

.00019 

.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


.00000 

.00000 

.00000 

.00000 


.18045 

.00004 

.00000 

.00000 


.18045 

.00004 

.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


.00000 

.00000 

.00000 

.00000 


.09022 

.00001 

.00000 

.00000 


.09022 

.00001 

.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00006 


1.00000 

1.00000 

1.00000 

1.00000 


.00000 

.00000 

.00000 

.00000 


.03609 

.00000 

.00000 

.00000 


.03609 

.00000 

.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1.00000 

1.00000 

1.00000 

1.00000 


1 .00000 
l.  00000 
1.00000 
1.00000 


.00000 

.00000 

.00000 

.00000 


.01203 

.00000 

.00000 

.00000 


.01203 

.00000 

.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 

1.00000  1.00000 


.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.13534 

.40601 

.67668 

.85712 

.94735 

.98344 

.99547 

.99890 

.99976 

.99995 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.40601 

.85712 

.98344 

.99890 

.99995 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.40601 

.45112 

.12631 

.01547 

.00105 

.00005 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN  FOR  BO 
.54416  .75458  .72114 

, TOTMU, & 
.00000 

AWP 

.00000 

.00000 

.00000 

2.00000 

2.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

PART  NO. 

4  *  *  * 

« 

* 

« 

* 

*  *  *  *  * 

*  *  *  * 

***** 

*  *  *  * 

SA  1  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  2  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  3  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  4  BO  CDF 

1.00000  1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

SA  5  BO  CDF 


1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  CDF 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

AWP  PF 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

AMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  PF 

.36788 

.36788 

.18394 

.06131 

.01533 

.00307 

.00051 

.00007 

.00001 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

TOTMU  CDF 

.36788 

.73576 

.91970 

.98101 

.99634 

.99941 

.99992 

.99999 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

l.  00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  CDF 

.36788 

.91970 

.99634 

.99992 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

1.00000 

BO  PF 

.36788 

.55182 

.07664 

.00358 

.00008 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

VAR/MEAN/VAR  TO  MEAN  FOR  BO 
.37902  .71617  .52924 

,  TOTMU, & 
.00000 

AWP 

.oooOo 

.00000 

.00000 

1.00000 

1.00000 

1.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 

.00000 
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Appendix  E 


Type  _I  Error  Tests 


Cases  with  N o  Type  1_  Error 
(HA  stock-0  and  HA  QPHA-1  ) 

The  first  five  cases  all  refer  to 


the  hierarchy  below; 


Case  1: 


Part 

HA 

SAs 


PART 

DATA 

MU 

Stock 

QPH  A 

1 

0 

1 

1 

0 

1 

HA  BO  RESULTS 

Mean  Variance 


EXACT  2.524  2.034 
APPROX  2.524  2.524 
%  error  0  24.1 


Case  2  : 


PART  DATA 


Part 

MU 

Stock 

QPHA 

HA 

1 

0 

1 

SAs 

1 

0 

2 

HA  BO  RESULTS 

Mean  Variance 

2.026  1.312 

2.026  2.026 
0  54.4 


EXACT 
APPROX 
X  error 


Part 

HA 

SAs 


PART  DATA 


MU 

Stock 

QPH  A 

1 

0 

1 

1 

1 

2 

HA 

BO  RESULTS 

Mean 

Variance 

EXACT 

1.497 

1 .330 

APPROX 

1.497 

1 .497 

%  error 

0 

12.6 

PART  DATA 


Part 

MU 

Stock 

QPHA 

HA 

2 

0 

1 

SAs 

3 

2 

2 

HA 

BO  RESULTS 

Mean 

Variance 

EXACT 

3.230 

2.701 

APPROX 

3.230 

3.230 

X  error 

0 

19.6 

Part 

PART 

MU 

DATA 

Stock 

QPHA 

HA 

4 

0 

1 

SA  1 

1 

0 

1 

SA  2 

7 

0 

1 

HA 

BO  RESULTS 

Mean 

Variance 

EXACT 

10.905 

10.377 

APPROX 

10.904 

10.444 

%  error 

"0 

.6 
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v  v  v  v  , 


Cases  six  through  nine  use  this  hierarchy: 


Case  6  *, 


PART 

DATA 

Part 

MU 

Stock 

OPHA 

HA 

1 

0 

1 

SAs 

1 

1 

2 

HA  BO  RESULTS 


Mean 

Variance 

EXACT 

1.879 

1.301 

APPROX 

1.879 

1.879 

X  error 

0 

44.4 

Case  7  : 

Part 

PART 

MU 

DATA 

Stock 

QPHA 

HA 

3 

0 

1 

SAs 

1 

2 

3 

HA 

BO  RESULTS 

Mean 

Variance 

EXACT 

3.345 

3.232 

APPROX 

3.345 

3.345 

Case  8  : 


PART  DATA 
Part  MU  Stock 

HA  2  0 

SAs  2  0 


QPH  A 
1 
1 


HA  BO  RESULTS 

Mean  Variance 


EXACT  5.731  3.513 
APPROX  5.731  5.731 
%  error  0  63.1 


Case  9: 


Part 
HA 
SA  1 
SA  2 
SA  3 
SA  4 
SA  5 


PART  DATA 
MU  Stock  QPHA 
2  0  1 

1  0  1 

2  0  1 

3  0  1 

4  0  1 

5  0  1 


HA  BO  RESULTS 


EXACT 
APPROX 
%  error 


Mean  Variance 

7.976  5.368 

7.976  7.976 

-0  48.6 


The  mean  errors  that  are  not  exactly  0  are  caused  by  the 
probability  distributions  being  truncated  in  the  computer 
programs . 


Cases  with  Type  1_  Error 

(HA  stock  not  equal  to  0  or  HA  QPHA  not  equal  to  1) 
Cases  10  through  13  use  the  following  hierarchy: 


Case  10: 

PART  DATA 

Part  MU  Stock  QPHA 

HA  1  2  1 

SAs  1  1  2 


HA  BO  RESULTS 

Mean  Variance 
EXACT  .2498  .3696 

APPROX  .2798  .4419 

Z  error  12.0  19.6 


Case  11: 

PART  DATA 
MU  Stock  QPHA 
1  0  3 

1  0  2 


HA  BO  RESULTS 

Mean  Variance 
EXACT  1.054  .1533 

APPROX  1.021  .2939 


Part 

HA 

SAs 


Case  12 


PART  DATA 
Part  MU  Stock 

HA  1  2 

SAs  1  0 


QPHA 

1 

2 


HA  BO  RESULTS 

Mean  Variance 


EXACT  .4370  .6109 
APPROX  .5571  .9223 
%  error  27.5  51.0 


Case  13: 


Part 

HA 

SA  1 
SA  2 


PART  DATA 
MU  Stock  QPHA 
3  5  1 

1  0  1 

5  0  1 


HA  BO  RESULTS 

Mean  Variance 
EXACT  3.175  6.655 

APPROX  3.188  6.748 

1  error  .4  1.4 


Cases  14  and  15  use  this  hierarchy 


Case  14: 


Part 

HA 

SAs 


PART  DATA 
MU  Stock  QPHA 
1  2  1 

1  1  2 


HA  BO  RESULTS 


EXACT 
APPROX 
X  error 


Mean  Variance 

.2998  .4364 

.3556  .5722 

18.6  31.1 


Case  15: 


Part 
HA 
SA  1 
SA  2 
SA  3 


PART 

DATA 

MU 

Stock 

QPHA 

1 

2 

1 

1 

0 

1 

2 

1 

2 

1 

0 

2 

HA  BO  RESULTS 

Mean  Variance 
EXACT  .7422  1.059 

APPROX  .8612  1.447 

%  error  16.0  36.6 
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PART 

DATA 

Part 

MU 

Stock 

QPHA 

HA 

1 

2 

1 

S  As 

1 

1 

2 

HA 

BO  RESULTS 

Mean 

Variance 

EXACT 

.  3395 

.4870 

APPROX 

.4191 

.6823 

X  error 

23.4 

40. 1 

Ca  se  17: 


PART 

DATA 

Part 

MU 

Stock 

QPHA 

HA 

1 

2 

1 

SAs 

1 

1 

2 

HA  BO  RESULTS 

Mean  Variance 
EXACT  .3717  .5264 

APPROX  .4714  .7732 

X  error  26.8  46.9 
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Appendix  F 


Type  2  Error  Tests 


Cases  with  N o  Type  2_  Error 

(Only  1  SA  per  HA  above  bottom  level  &  no  type  1  error) 


Case  1 


PART  DATA 


Part 

1 

2 

3 

A 


Stock 

0 

0 

2 

3 


QPHA 

1 

1 

4 

2 


EXACT 
APPROX 
%  error 


RT  1  BO  RESULTS 

Mean  Variance 
5.688  5.319 

5.688  5.688 

r  0  6.9 


Case  2: 


Part 

1 

2 

3 

4 

5 


PART  DATA 
MU  Stock 


QPHA 

1 

1 

1 

1 

1 


PART 

1  BO  RESULTS 

Mean 

Variance 

EXACT 

5.448 

4 .508 

APPROX 

5 .448 

5.448 

X  error 

0 

20.9 

PART  DATA 

Part  MU  Stock  QPHA 
110  1 

2  3  0  1 

3  10  1 

4  2  3  2 

5  10  3 

6  3  2  4 


PART  1  BO  RESULTS 

Mean  Variance 

5.936  5.201 

5.936  5.936 

0  14.1 


EXACT 
APPROX 
X  error 


Cases  with  Type  £  Error 

Cases  four  through  six  use  the  following  hierarchy 


Case  4  : 

PART  DATA 

Part  MU  Stock  QPH A 
all  2  0  1 


PART  1  BO  RESULTS 

Mean  Variance 
EXACT  7.856  5.144 

APPROX  7.985  7.927 

I  error  1.6  54.1 


Case  5  : 


Part 

PART 

MU 

DATA 

Stock 

QPHA 

1 

4 

0 

1 

1  St 

2 

0 

1 

PART 

1  BO 

RESULTS 

Mean  Variance 


EXACT  9.856  7.145 
APPROX  9.987  9.987 
I!  error  1.3  39.8 


Part 

1-3 

4-7 


PART  DATA 
MU  Stock  QPH A 
.5  0  1 

.5  1  1 


PART  1  BO  RESULTS 

Me  an  Variance 


EXACT  1.634  1.347 
APPROX  1.628  1.628 
Z  error  -.4  20.9 


The  next  case  uses  this  hierarchy: 


Case  7: 


Part 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 


PART 

DATA 

MU 

Stock 

QP11A 

2 

0 

1 

1 

0 

1 

3 

0 

1 

1 

0 

1 

1 

0 

1 

2 

0 

1 

1 

0 

1 

2 

0 

1 

3 

0 

1 

2 

0 

1 

A/C  BO  RESULTS 

Mean  Variance 
7.587  3.942 

7.782  5.138 


EXACT 

APPROX 


Appendix  G 

Sensitivity  Analysis  Results 


Effect  of  Increasing  Depth  of  Indenture 


For  the  single-level  hierarchy: 


PART  DATA 

Part  MU  Stock  QPHA 

12  3  2 

2  1.5  1  2 


A/C  BO  RESULTS 

Mean  Variance 
EXACT  .6088  .4239 

APPROX  .6088  .4239 

%  error  0  0 


For  the  ''two-level”  hierarchy: 

A/C 


odd  #s 
even  #s 


PART  DATA 


MU 

Stock 

QP 

2 

3 

2 

1 . 5 

1 

2 

rwwv.*' 


A/C  BO  RESULTS 

Mean  Variance 
EXACT  .9422  .5245 

APPROX  .9557  .5683 

U  error  1.4  8.4 


For  the  "three-level"  hierarchy: 


Repeating  the  same  sequence  for  a  different  data  set  gives 
the  following  results: 

For  the  " s i ng le- 1  eve  1 “  hierarchy: 


For  the  "two -lev el"  hierarchy: 


PART  DATA 

Part  MU  Stock  QPHA 
all  2  3  1 


A/C  BO  RESULTS 

Mean  Variance 
EXACT  .4126  .6568 

APPROX  .4126  .6568 

1  error  0  0 


PART  DATA 
MU  Stock  QPHA 
2  3  1 


-  *  .*  V  <*  v  ’  •  V  V  *.  .  ,  V 


PART  DATA 

Part  MU  Stock  QPH A 
all  2  3  1 


A/C  BO  RESULTS 

Mean  Varlan 
EXACT  1.405  2.643 

APPROX  1.091  1.685 

%  error  -22.3  -36.2 


PART  DATA 

Part  HU  Stock  QPH A 

all  1.5  2  2 


EXACT 
APPROX 
%  error 


A/C  BO  RESULTS 

Mean  Variance 
1.053  .4226 

1.106  .4822 

r  5.0  14.1 


For  an  SA/HA  of  five: 


Part 

all 


PART  DATA 
MU  Stock  QPH A 

2  3  1 


A/C 


EXACT 
APPROX 
X  error 


BO  RESULTS 
Mean  Variance 

1.654  2.036 

1.562  1.794 

-5.6  -11.9 


f  ive : 


PART  DATA 

Part  MU  Stock  QPH A 
all  2  3  1 


A/C  BO  RESULTS 

Mean  Variance 
1.876  2.056 

1.804  1.877 

r  -3.8  -8.7 


EXACT 

APPROX 


A/C 


PART  DATA 

Part  MU 

Stock 

QPHA 

all  1 

0 

1 

A/C 

BO  RESULTS 

Mean 

Variance 

EXACT 

A  .533 

1.458 

APPROX 

4.924 

2.228 

%  error 

8.6 

52.8 

For  an  SA/HA  of  five: 

PART  DATA 

Part  MU 

all  1 

Stock 

0 

QPHA 

1 

A/C 

BO  RESULTS 

EXACT 

APPROX 

7.  error 

Mean 

7.231 

7.992 

10.5 

Variance 

2.743 

5.229 

90.6 

145 

L_*V‘  " 


PART 

DATA 

Part 

MU 

Stock 

QPHA 

all 

1 

3 

1 

A/C  BO  RESULTS 

Mean 

Variance 

EXACT  .1354 

.1940 

APPROX  .1219 

.1678 

Z  error  -10.0 

-13.5 

five: 


PART  DATA 

Part  MU  Stock  QPH A 

all  1  3  1 


A/C  BO  RESULTS 

Mean  Variance 
EXACT  .0632  .0940 

APPROX  .0556  .0787 


Effect  o f  Stock  Level  and  QPHA 

This  hierarchy  with  all  MUs  equal  to  two  was  used  to  gen¬ 
erate  the  results  in  the  following  table: 


Stock 

Level 

QPHA 

mean 

1 

var . 

mean 

2 

var . 

mean 

3 

var. 

APPROX 

9.569 

6.289 

3.145 

1.099 

1 .941 

.4306 

0 

EXACT 

9.119 

4.131 

2.946 

.6146 

1.848 

.3105 

X  error 

4.9 

52.2 

6.8 

78.8 

5.0 

38.7 

APPROX 

6.222 

4.733 

2 .195 

.8614 

1.433 

.  3792 

1 

EXACT 

6.122 

4.122 

2.071 

.6156 

1.353 

.2840 

X  error 

1.6 

14.8 

6.0 

40.0 

5.9 

33.5 

APPROX 

3.061 

3.234 

1.248 

'  1 

.9184 

.3542 

2 

EXACT 

3.215 

3.825 

1.201 

.8900 

.3027 

X  error 

-A. 8 

-15.5 

3.9 

1 

3.2 

17.0 

APPROX 

.9043 

1.420 

.4271 

.4061 

.2860 

3 

EXACT 

1.066 

1.948 

.4183 

.3952 

.2772 

mi 

-27.1 

2.1 

2.8 

3.2 

APPROX 

.2060 

.3516 

.1406 

.1496 

.1282 

mi 

4 

EXACT 

.2465 

.4648 

.1418 

.151 1 

.1278 

X  error 

-16.4 

-24.4 

-.9 

-1.0 

.  3 

mil 

APPROX 

.0500 

.0829 

.0382 

.0421 

m 

.0358 

5 

EXACT 

.0548 

.0943 

.0385 

.0425 

.0358 

X  error 

-8.7 

-12.1 

-.7 

-.8 

.06 

.03 

This  hierarchy  with  all  MUs  equal  to  2.7  was 
erate  the  results  in  the  following  table: 


QPHA 


2 


St  oc 
Leve 


APPROX 

0  EXACT 

X  error 

7.185  5.137 
7.085  3.906 
1.4  31.5 

APPROX 

1  EXACT 

X  error 

5.420  3.872 
5.355  3.361 
1.2  15.2 

APPROX 

2  EXACT 

X  error 

3.779  2.842 
3.774  2.812 
.1  1.1 

APPROX 

3  EXACT 

%  error 

2.333  2.168 
2.369  2.319 
-1.5  -6.5 

APPROX 

4  EXACT 

X  error 

1.197  1.539 
1.231  1.650 
-2.8  -6.7 

3.193 

3.124 

2.2 


2.551 

2.499 


2 


1.926 

1.859 

3.6 


1.328 

.317 

8 


.7704 
.  7700 
.  1 


.8176 

.6357 

28.6 


7039 

5934 

8.6 


.6268 
.  5735 
9.3 


.5633 
.  5457 
3.2 


.5108 

.5088 


mean 

3 

var  . 

2 .191 

.  3253 

2.152 

.2675 

1.8 

21.6 

1.795 

.3507 

1.765 

.3195 

1.7 

9.8 

.3076 

1.378 

.2869 

7.2 

1.052 

.  2678 

1.043 

.  2574 

.9 

4.0 

.6631 

.3191 

.6602 

.3172 

.4 

.6 
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To  accomplish  this  two  Fortran  computer  programs  were 
developed,  one  to  compute  the  expected  number  of  not  mission 
capable  aircraft  using  the  accurate  mathematical  calculations, 
the  other  to  compute  the  same  figure  using  the  Poisson  approx¬ 
imation.  These  programs  were  used  to  evaluate  the  approxi¬ 
mation  caused  error  for  different  parts  hierarchies  and  data 
sets. 

The  analysis  identified  two  distinct  causes  for  the 
error  induced  by  the  approximation.  These  sources  of  error 
were  confirmed  by  running  test  cases  specifically  tailored  to 
eliminate  the  error-causing  characteristics  and  noting  that 
no  approximation  error  resulted. 

The  approximation  error  was  found  to  fluctuate  in  sign 
and  magnitude  for  different  cases.  Sensitivity  analysis  was 
performed  to  identify  the  sensitive  parameters. 
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